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STORM  SURGE  SIMULATION  IN  TRANSFORMER  COORDINATES 
Volume  I.  Theory  and  Application 
by 

John  J.  Wanstrath,  Robert  E.  Whitaker, 

Robert  0.  Reid,  and  Andrew  C.  Vastano 

I.  INTRODUCTION 

1 . Background. 

Storm  surges  are  transient  fluctuations  in  the  sea  level  induced 
by  atmospheric  disturbances , notably  those  due  to  extra-tropical 
storms  and  hurricanes  and  to  a less  frequent  extent  pressure  jumps 
associated  with  line  squalls.  The  rise  of  the  water  and  circulation 
caused  by  a hurricane  can  be  considerable  and  is  of  special  practi- 
cal importance  with  respect  to  loss  of  lives  and  property  not  only 
adjacent  to  the  coast  but  also  well  inland.  Statistical  studies  of 
hurricanes  of  record  provide  a means  of  predicting  the  surge  height 
along  an  open  coast.  The  empirical  formulas  developed  from  these 
studies  relate  the  maximum  expected  surge  height  to  meteorological 
parameters  and  effective  coastal  bathymetry  (Donn,  1958).  However, 
all  such  studies  do  not  provide  the  time  history  or  even  the  time 
scale  of  the  forcing  function  which  is  necessary  as  input  for  bay- 
response  studies  (Reid  and  Bodine,  1968). 

More  recently,  time-dependent  models  based  upon  the  physics  of 
the  storm  surge  phenomena  have  been  developed  to  study  the  generation 
or  modification  of  the  surge  as  it  leaves  deep  water  and  moves  over 
the  Continental  Slope  and  Shelf.  These  models,  like  the  one  proposed 
herein,  involve  the  vertically  integrated  equations  of  motion  and 
mass  continuity.  The  greatest  difficulty  in  utilizing  these  models 
has  been  the  manner  in  which  the  shoreline  has  been  portrayed  and  the 
application  of  realistic  boundary  conditions  at  the  specified  shore. 
.Jelesnianski  (1967,  1972)  takes  the  shoreline  as  a vertical  plane  of 
infinite  height,  thereby  facilitating  the  mathematical  representa- 
tion of  the  shore  boundary.  More  general  portrayal  of  the  shoreline 
is  achieved  by  the  schemes  of  Miyazaki  (1963)  and  Platzman  (1963)  in 
which  the  coastline  is  represented  as  a series  of  straight-line  seg- 
ments connected  at  right  angles.  Specifying  the  shoreline  in  this 
stairstep  manner  results  in  greater  numerical  programing  complexity 
for  the  shore  boundary  condition  in  that  the  algorithm  must  possess 
the  ability  to  search  for  and  substantiate  the  location  of  land.  A 
more  serious  objection  is  that  this  approximation  may  inject  spurious 
oscillations  into  the  calculations.  This  adverse  effect  is  more  than 
academic  since  the  concentration  of  energy  is  near  its  maximum  at  the 
coast,  and  further,  this  is  precisely  the  region  from  which  water 
level  observations  are  usually  available  for  model  verification. 


II 


Clearly,  these  difficulties  are  avoided  or  minimized  if  a curvi- 
linear coordinate  system  is  employed  which  will  map  the  interior 
region  bounded  by  the  actual  coastline,  the  seaward  boundary  (taken 
as  the  180-meter  depth  contour)  and  two  lateral  straight  lines.  Thus, 
the  transformation  is  such  that  the  coast,  as  well  as  the  180-meter 
isobath,  is  mapped  as  a straight  line  in  the  image  plane.  In  order 
to  minimize  the  associated  modification  of  the  differential  equations 
governing  the  storm  surge,  an  orthogonal  curvilinear  system  is 
employed  in  the  present  study,  the  desired  mapping  being  carried  out 
via  a conformal  transformation. 

2 . Objectives. 

The  objectives  are  twofold:  (a)  to  develop  a method  which  will 

conformally  map  the  interior  region  bounded  by  two  arbitrarily  shaped 
curves  each  beginning  and  ending  on  parallel  straight  lines  drawn 
between  the  two  curves;  (b)  to  develop  and  verify  a numerical  model 
of  a hurricane-induced  surge  on  the  Continental  Shelf  employing  the 
transformed  equations  of  motion  in  the  orthogonal  curvilinear 
coordinate  system. 


II. 


CONFORMAL  MAPPING 


1 . Development. 

The  selection  of  appropriate  conformal  mapping  and  error  functions 
is  of  principal  concern  in  achieving  the  first  stated  objective.  We 
desire  to  map  a region  R of  the  Z-plane  into  a rectangle  in  the  c- 
plane  (Figure  1),  in  which  the  seaward  boundary  and  coastline  curves 
are  specifically  transformed  into  the  image  plane  as  constant  values 
of  n • Furthermore,  the  curves  in  the  negative  x region  are  re- 
quired to  be  the  mirroT  image  of  those  in  the  positive  x region. 

As  an  artifice  to  assure  that  the  lateral  boundaries  of  the  mapped 
region  represent  straight  parallel  lines  normal  to  the  shoreline  and 
the  seaward  boundary,  symmetry  about  one  of  these  boundaries  (x  = 0) 
is  imposed  and  the  whole  range  in  x (-A  to  A)  is  considered  to  be  one 
wavelength  of  a periodic  function.  Thus  only  the  range  0 < x < A 
in  the  top  panel  of  Figure  1 corresponds  to  the  real  region  of  the 
shelf. 


The  conformal  mapping  relation  is  taken  in  the  form: 

Z = F(0  , 

where 


(1) 


Z = x + iy  , (2) 

C * S ♦ in  , (3) 

and 

i « • (4) 

An  appropriate  form  of  the  transformation  for  the  mapping  considered 
in  Figure  1 is: 

N 

F(C)  = PQ  + QqG  + l (Pn  cos  nkc  + sin  nke)  , (5) 

n=l 


where 

k = tt/A  , (6) 

and  A is  half  the  horizontal  extent  of  the  region  in  the  Z , or  c 
plane.  The  coefficients  PR  and  Qn  , n=0,l...N,  are,  in  general, 
complex  and  independent  of  t or  n . Let 


A + 
n 


iB 


(7) 
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Figure  1.  Conformal  mapping  planes  for  transforming  the  interior  region  of 
the  Z-plane  into  the  interior  region  of  the  c-plane. 


r 

and 

% - Cn  * i0n  • (»> 

where  the  coefficients  A„  , Bn  , Cn  , and  Dn  are  real  quantities. 
However,  the  symmetry  condition  about  x = 0 (or  £ = 0)  requires  that 
all  Ajj  and  Dn  be  zero.  The  resulting  relations  for  x and  y in 
terms  of  £ and  n are: 

N 

x(£,o)  = Cq£  + £ (B^  sinh  nkn  ♦ Cn  cosh  nkq)  sin  nk£  , (9) 

n=l 

and 

N 

y(£»n)  = Bq  + CQn  + £ (B^  cosh  nkn  + sinh  nkn)  cos  nk£  . (10) 

n=l 

The  condition  that  the  range  of  x and  £ be  the  same  requires  that 
CQ  = 1.  The  remaining  N values  of  Cj,  and  N+l  values  of  Bj,  are 
determined  by  matching  the  coastal  and  seaward  boundary  curves  at 
n = 4 6,  respectively,  3 also  being  a parameter  to  be  determined. 

The  coordinates  X,Y  of  the  given  coast  or  seaward  boundary 
curves  are  specified  parametrically  in  terms  of  arc  length  measured 
along  each  curve  from  some  fixed  point.  The  functions  Xs  , Ys  , 

Xc  , and  Yc  are  single-valued  functions  of  this  parameter  where  the 
superscript  s or  c represents  the  seaward  boundary  or  coastline 
curve,  respectively.  This  property  is  essential  since  a Fourier 
series-type  representation  is  employed  in  determining  the  coefficients. 
The  problem  is  to  determine  3 , B0  , and  the  set  of  coefficients, 

Bn  and  Cn  , for  a given  N such  that  equations  (9)  and  (10)  give  a 
best  fit  to  the  given  curves,  in  the  sense  of  minimizing  an  appropriate 
mean  square  error  function.  Since  the  specified  curves  are  not  known 
directly  in  terms  of  £ , but,  rather  in  terms  of  arc  length,  the 
bicurve  fitting  equations  require  an  iterative  process  starting  from 
some  initial  estimate  of  arc  length  in  terms  of  £ for  each  curve. 

It  is  found  that  a convergent  iterative  procedure  results  if  one 


chooses  the  following  as  the  i+1  approximation  of  the  coefficients: 


(13) 


b = fy  a - r * )/ (t  « - ($  )2)  , 

n n n n n'  n n n'  1 


C„  - <r„  T„  - rn  V'<T»  “n  ' <V!>  • ',4> 


where  i-1,2...  denotes  the  iteration  number,  A is  the  transform 
generated  arc  length  and 


Tn  = ~ — 2~^ — sinh2  nkB*+1  + — cosh2  nkBi+1  , (15) 


(N^'-W**1)  (Wy^-W^’1)  .+1 

4>  = 2- t-  * sinh  2nk6  , 

n 4 


wc,i+ws,i« 

— sinh2  nkB*+1  + — — ■— — cosh2  nkBi  + 1 ,(17) 


Tn  = S1-— yk-S [W^’1  / (Xc (A1)  - C)  sin  nk£  dC 


. A 

- W*’1  / (Xs (A1)  - a sin  nkc  d£  ] 
* o 


, . „i+l  . A 

cosh  nkB  ru,c,i  , „c,.i.. 

+ r — [W  ’ / Y (A  ) cos  nk£  d£ 

A / o 


. A 

♦ W*’1  / YS (A1)  cos  nkC  dC]  , 
y o 


X 


(24) 


Ey’i+1  = lx  - y1+1(^.+B)]2  d*  * 

-X 

The  transform  generated  arc  length  and  weighting  factors  are  deter- 
mined relative  to  the  results  of  the  previous  iteration  by: 

WS  = ES/E  , (25) 

x x 

W * = E^/E  , (26) 

Wx  = Ex/f  ' (2?) 

Wy  = E^/E  , (28) 

where  E is  the  average  value  of  equations  (21)  through  (24).  The 
arc  length  functions  are  given  by: 


A(£+,-6)  = yS  J [(ff)2  + (ff)2]15  d?  , (29) 

o 

and 


A(C+,8)  = YC  f [(||)2  ♦ (f^)2]15  dC  , (30) 

o 

for  0 < £+  < X and  yS  and  are  prorating  factors  taken  such 

that  Alx,n)“  is  exactly  the  known  length  of  the  given  curve.  Two 
levels  of  approximation  are  involved  in  the  error  function.  The  ith 
approximation  of  these  quantities  is  involved  in  the  interpolation  of 
the  given  curves.  This  approximation  is  mandatory  in  deriving  a 
deterministic  form  for  these  quantities  relative  to  the  i+1  level. 

The  iterative  procedure  is  initiated  by  setting  all  coefficients 
to  zero,  all  weighting  factors  to  unity,  and  equations  (29)  and  (30) 
are  given  by  ysS+  and  YcC^  » respectively.  The  first 
estimates  of  6 , Bq  and  the  Fourier-type  coefficients 
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arc  determined  from  equations  (11)  through  (14),  respectively.  Hav- 
ing the  first  approximation  for  these  quantities,  the  transform- 
generated  arc  length  for  each  curve  and  weight  factors  may  be  deter- 
mined thus  permitting  a second  approximation.  This  procedure  is 
repeated  for  given  N until  6 , BQ  and  the  transform  coefficients 
have  converged  as  indicated  by  successive  values  of  the  error  func- 
tion. The  error  function  approaches  a constant  value  governed  by  N 
and  the  inherent  errors  in  estimating  the  various  integrals  by 
numerical  methods.  The  trapezoidal  rule  employing  M values  of  £ 
at  equal  increments  over  the  interval  0 < £ £ x is  used  to  estimate 
all  integrals.  Given  the  value  of  £ , the  appropriate  values  of 
X , Y , Xc  and  Yc  are  computed  by  an  interpolational  subroutine 
called  CURV,  which  was  obtained  from  Dr.  Alan  Cline,  National  Center 
for  Atmospheric  Research.  This  routine  is  based  on  a numerical 
analogy  to  a spline  under  tension  (for  example,  see  Schweikert, 

1966).  In  all  applications  of  the  conformal  mapping  procedure, 
the  total  number  of’  ? erms  of  the  series,  N , is  taken  as  M or 
smaller.  The  limitation  on  N is  imposed  since,  for  a determinate 
system,  the  total  number  of  Fourier-type  coefficients  cannot  exceed 
the  number  of  data  points.  The  selection  of  N is  governed  by  the 
desired  degree  of  fit  of  the  functions  Xs  , Ys  , Xc  and  Y0 
The  procedure  outlined  above  follows  the  rudiments  of  a conformal 
mapping  technique  developed  by  Reid  and  Vastano  (1966) . 

2.  A££l  ications . 


Three  regions  of  the  Continental  Shelf  of  the  Gulf  of  Mexico  and 
two  regions  of  the  eastern  seaboard  of  the  United  States  have  been 
mapped.  These  regions  are: 


(a)  western  gulf  coast: 

(b)  central  gulf  coast: 

(c)  eastern  gulf  coast: 

(d)  lower  east  coast: 

(e)  upper  east  coast: 


Laguna  Madre,  Mexico,  to  Marsh 
Island,  Louisiana; 

Matagorda  Bay,  Texas,  to  Timbal ier 
Bay,  Louisiana; 

Atchafalaya  Bay,  Louisiana,  to 
Apalachee  Bay,  Florida; 

Cape  Kennedy,  Florida,  to  Pamlico 
Sound,  North  Carolina; 

Pamlico  Sound,  North  Carolina,  to 
Penobscot  Bay,  Maine; 


and  are  shown  in  Figures  2 through  6.  The  top  and  bottom  curve  in 
each  figure  represents  the  coastline  and  the  180-meter  isobath, 
respectively.  The  latter  is  taken  to  correspond  approximately  to 
the  shelf  break,  seaward  of  which  the  depths  increase  abruptly;  the 
180-meter  contour  may  not  always  be  the  most  appropriate  one  for  the 
outer  limit  of  the  Continental  Shelf.  The  chart  scale,  orientation 
of  the  region,  and  location  of  cities  and  bays  along  the  coast  are 
indicated  on  each  figure.  The  x's  represent  the  location  where 


Figure  3.  Coastline  (.top)  and  seaward  boundary  curve  (bottom),  central  gulf  coast 


Coastline  (top)  and  seaward  boundary  curve  (bottom),  eastern  gulf  coast 


Figure  5.  Coastline  (.top)  and  seaward  boundary  curve  (bottom),  lower  east  coast. 


the  coast  and  the  180-meter  contour  were  digitized  from  an  overlay 
of  bathymetric  charts  for  the  regions  concerned.  The  smooth  curve  is 
the  result  of  applying  the  numerical  spline  interpolation  routine  to 
these  points.  In  reading  the  discrete  positions,  both  the  coastline 
and  the  180-meter  isobath  have  been  smoothed  subjectively  to  suppress 
variations  with  scale  lengths  less  than  the  grid  scale  ultimately 
employed  in  the  surge  calculations.  For  example,  narrow  entrances 
to  bays  are  replaced  by  a fictitious  coastline  across  the  entrance, 
and  the  cusp-shaped  features  are  smoothed. 

A method  was  adopted  which  minimizes  the  amount  of  computer  time 
required  to  solve  equations  (9)  and  (10)  by  hastening  the  convergence 
of  the  iterative  procedure.  The  method  consisted  of  increasing  N 
at  selected  iteration  intervals.  The  iterative  procedure  used  to 
determine  the  transformation  coefficients  was  terminated  if  the  con- 
vergence criterion  was  satisfied  or  if  the  available  computer  storage 
was  exceeded  because  of  increasing  N . The  convergence  criterion 
was  a mean  variance  of  less  than  1 square  kilometer  between  the 
transform-generated  curves  and  that  specified.  However,  for  develop- 
ment purposes,  the  iterative  procedure  for  the  western  and  central 
gulf  coast  and  lower  east  coast  areas  was  continued  beyond  this 
criterion  to  obtain  a better  fit. 

It  is  convenient  to  discuss  the  mapping  of  the  western  and  central 
gulf  coast  regions  at  the  end  of  this  section  since  additional  testing 
of  the  solution  to  the  mapping  equations  was  performed  with  these 
areas.  Figure  4 shows  the  shelf  region  for  the  eastern  gulf  coast 
where  the  Mississippi  Delta  is  shown  as  the  shaded  area.  Considera- 
tion of  the  numerical  time  step  for  the  surge  algorithm  was  the  prin- 
ciple reason  why  the  coastline  to  be  mapped  did  not  follow  the  delta. 
The  reduction  in  the  time  step  if  the  actual  coastline  had  been 
followed  would,  probably,  have  been  at  least  tenfold  with  respect  to 
the  one  used  in  the  simulation  of  the  storm  surge  induced  by  Hurri- 
cane Camille  (Section  V).  Since  the  delta,  or  really  the  levee 
adjacent  to  the  Mississippi  River,  has  a controlling  influence  on 
the  circulation  and  surge  caused  by  the  hurricane,  this  geographical 
feature  was  included  in  the  surge  model  as  a wall  protruding  from  the 
coast.  The  fit  of  the  transform-generated  coast  and  seaward  boundary 
curves  with  respect  to  that  specified  after  1,  20,  40,  80,  and  160 
iterations  is  shown  in  Figures  7 through  11,  respectively.  The  suc- 
cessively better  agreement  of  the  mapped  curves  with  respect  to  that 
specified  is  obvious  from  the  convergence  table  for  the  gulf  coast 
region  of  Atchafalaya  Bay  to  Apalachee  Bay  (Table  1).  The  table 
shows  the  variance  and  convergence  behavior  of  selected  transforma- 
tion coefficients  as  the  iteration  and  N increases.  The  coast 
variance  in  X , the  third  column  of  Table  1,  is  calculated  by: 

. M 

jCPT  l |XC(&  Uj'B))  - x^.B)!2  (31) 


25 


curve  after  one  iteration 


Figure  8.  Transform-generated  coastline  and  seaward  boundary  curve  after  20  iterations 
(solid)  and  that  specified  (dashed),  eastern  gulf  coast. 


eaward  boundary  < 
, eastern  gulf  c 


Table  1.  Convergence  table  for  the  gulf  coast  region  of  Atchafalaya  Bay  to  Apalachee  Bay 


Similar  expressions  are  employed  for  the  other  variances.  Appendix 
A contains  the  transformation  coefficients  used  to  produce  Figure  11 

The  successful  application  of  the  mapping  equations  to  the  other 
shelf  regions  was  accomplished.  Tables  2 through  5 clearly  indicate 
the  successive  convergence  of  the  mapped  curves  with  respect  to  that 
specified  for  the  lower  and  upper  east  coast  and  the  western  and 
central  gulf  coast,  respectively.  Figures  12  through  15  show  the 
fit  of  the  transform- generated  curves  and  that  specified  after  the 
last  iteration  for  the  above  regions.  Appendixes  B and  C contain 
the  mapping  coefficients  used  to  produce  Figures  12  and  13,  respec- 
tively. The  bottom  part  of  Tables  4 and  5 show  the  results  of  addi- 
tional testing  of  the  conformal  mapping  equations  (9)  and  (10)  with 
the  western  and  central  gulf  coast  regions.  Another,  less  general, 
solution  of  equations  (9)  and  (10)  is  possible  if  one  minimizes  the 
least  square  error  function  defined  only  in  terms  of  the  Y inte- 
grals. This  solution  for  the  transformation  coefficients,  hereafter 
referred  to  as  an  alternate  solution,  may  be  obtained  from  equations 

(13)  through  (19)  with  W*  = Wx  = ° ' 7,10  testing  procedure  was  to 

continue  the  iterative  process  as  outlined  in  the  previous  section 
with  the  initial  approximation  for  the  coefficients  being  those 
values  determined  from  the  80th  iteration.  The  alternate  solution 
as  applied  to  the  western  and  central  gulf  coast  regions  was  stable 
and,  moreover,  provided  a better  fit  with  N = 110  than  the  more 
general  one.  The  fit  of  the  mapped  curves  with  respect  to  that 
specified  after  1 and  40  (or,  81  and  120)  iterations  with  the  alter- 
nate solution  is  shown  in  Figures  16  and  17,  respectively,  for  the 
western  gulf  coast  and  Figures  18  and  19,  respect i vely,  for  the  cen- 
tral gulf  coast.  Appendixes  I)  and  E contain  the  transformation  co- 
efficients used  to  produce  Figures  17  and  19,  respectively.  In  test 
ing  with  the  other  three  shelf  regions,  the  alternate  solution  was 
nonconvergent  in  that  the  successive  values  of  the  error  function 
do  not  decrease  or  approach  a constant  as  outlined  in  the  previous 
section.  The  reason  for  this  behavior  has  not  been  investigated. 

The  total  variance  of  x and  y from  a linear  transform  of  the 
curvilinear  coordinates  may  be  defined  as: 


where 

e 

2 1 
°x  " 26A 

'6 

and 


x y 
X 

, [x  - S]2  d£  dn  ; (33) 

o 


32 


Table  2.  Convergence  table  for  the  east  coast  region  of  Camp  Kennedy  to  Pamlico  Sound. 
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Table  3.  Convergence  table  for  the  east  coast  region  of  Pamlico  Sound  to  Penobscot  Bay 


Table  4.  Convergence  table  for  the  gulf  coast  region  of  Laguna  Madre  to  Marsh  Island. 
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lerated  coastline  and  seaward  boundary  curve  after  135  iterations 
hat  specified  (dashed),  upper  east  coast. 


Transform-generated  coastline  and  seaward  boundary  curve 
(solid)  and  that  specified  (dashed),  western  gulf  coast. 


Figure  15.  Transform-generated  coastline  and  seaward  boundary  curve  after  80  interations 
(solid)  and  that  specified  (dashed),  central  gulf  coast. 


Transform-generated  coastline  and  seaward  boundary  curve  after  one  iteration 
with  the  alternate  solution  (solid)  and  that  specified  (dashed),  western 
gulf  coast. 


ansf< 


msform-generated  coastline  and  seaward  boundary  curve  after  one 
:h  the  alternate  solution  (solid)  and  that  specified  (dashed) , ce 


jure  19.  Transform-generated  coastline  and  seaward  boundary  curve  after  40  iterations  with 
the  alternate  solution  (solid)  and  that  specified  (dashed),  central  gulf  coast. 
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[y  - (Bq  + n)  ]2  df,  dn  . (34) 

-8  o 

These  expressions  are  measures  of  the  "curvilinearity"  of  the  trans- 
forms x(£,n)  and  v(£,n)  . From  equations  (9)  and  (10),  it  can  be 
shown  that 

o2  = l on2  , (35) 

n=l 

where  the  contribution  to  the  total  curvilinearity  variance  due  to 
the  n^  term  is: 


(B  2 + C 2)  sinh  2nk8 
2 _ n n_ 

°n  4nk8 


(36) 


The  value  of  o 2/o2  is  a more  meaningful  measure  of  the  relative 
importance  of  the  various  terms  in  the  transformation  than  the  raw 
Bn  , C coefficients.  The  curvilinearity  variance  spectrum  given 
by  an2/o2  as  a function  of  n for  the  final  transform  of  the  five 
regions  is  shown  in  Figure  20.  The  values  of  a2  and  o/X  for' the 
five  regions  are  indicated  in  the  figure.  Each  of  the  curves  indi- 
cates a general  inverse  power  law  trend  of  about  4th  degree.  The 
variability  of  the  spectrum  in  all  but  the  middle  curve  is  more  pro- 
nounced in  the  higher  harmonics  than  in  the  lower. 

3.  [.imitations . 

Testing  of  the  conformal  mapping  equations  was  conducted  with 
coastlines  that  varied  in  the  degree  of  smoothness  and  regions  that 
varied  in  X . Another  version  of  the  coastline  from  Cape  Kennedy 
to  Pamlico  Sound  is  shown  in  Figure  21.  Comparing  Figure  5 with  this 
version,  it  can  be  seen  that  there  are  two  small  differences  in  the 
form  of  the  coastline.  Figure  21  shows  the  narrowing  of  the  Con- 
tinental Shelf  in  the  area  south  of  Cape  Kennedy.  Another  differ- 
ence is  the  slightly  sharper  vertex  at  Cape  Kennedy,  Cape  Fear,  and 
Cape  Lookout.  The  (general)  solution  of  the  mapping  equations  proved 
unstable  for  this  situation.  It  is  not  known  which  change  in  the 
coastline  permitted  the  convergence  of  the  iterative  procedure  for 
Figure  5. 


One  test  of  the  mapping  equations  was  conducted  on  the  region 
from  Matagorda  Bay  to  Apalachee  Bay  (Figure  22).  This  region  is  the 
result  of  joining  Figures  3 and  4 and  represents  in  this  study  the 
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Figure  21.  Coastline  (top)  and  seaward  boundary  curve  (bottom),  lower  east  coast. 


longest  horizontal  extent  of  a region  to  be  mapped.  The  mean  vari- 
ance between  the  mapped  curves  and  that  specified  after  120  itera- 
tions with  n=140  was  4.4  square  kilometers.  The  level  of  vari- 
ance achieved  indicates  that  the  bicurve  fitting  equations  are 
marginally  successful  in  mapping  this  enlarged  region.  The  change 
in  the  variance  level  over  the  last  10  iterations  and  limits  on 
computer  time  were  the  bases  for  terminating  the  iterative  procedure. 
The  convergence  criteria  can  be  satisfied  either  by  a shorter  hori- 
zontal extent  or  smoother  curves  for  a longer  region. 
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III.  THF.  STORM  SURGE  EQUATIONS  IN  THE 
SHELF  COORDINATE  SYSTEM 


1 . Stretched  Shelf  Coordinate  System. 

Consider  the  transform-generated  coastline  and  seaward  boundary 
curves  as  shown  in  Figure  23.  The  orthogonal  curvilinear  mesh  asso 
ciated  with  the  coordinates  (£,n)  is  designated  the  shelf  coordi- 
nate system.  Although  the  surge  calculations  can  be  performed  in 
the  shelf  coordinate  system  (Figure  23b),  a problem  remains  in 
obtaining  the  desired  spatial  resolution  with  the  fewest  possible 
computational  points.  The  computing  grid  employed  (Figure  23c)  is 
the  result  of  a second  transformation  which  preserves  the  orthogo- 
nality property  but  allows  independent  stretching  of  £ and  of  n 
The  grid  resulting  from  the  second  transformation  is  termed  the 
stretched  shelf  coordinate  system  (S*,T*). 

The  stretched  shelf  coordinate  system  is  generated  by  indepen- 
dently transforming  the  £ and  n axes  in  the  following  manner: 

(a)  Given  the  hurricane  storm  track,  the  location  where 
the  storm  crosses  the  coast  and  the  coastline  con- 
figuration, the  values  of  £ along  the  coastline 
?re  determined  which  will  produce  a constant  rela- 
tively fine  increment  of  coastline  arc  length,  Sp  . 

In  this  area  of  prime  interest,  line  BC  in  Figure^ 

23(a),  the  constant  increment  of  arc  length  is 
equal  to  AS*  . However,  regions  AB  and  CD  show 
that  for  the  same  AS*  as  above,  there  is  a rela- 
tive expansion  of  the  increment  of  the  coastline 
arc  length.  The  functional  relationship  between  £ 
and  S*  is: 

S*  = S* (Sp (£) ) , (37) 

where  the  expansion  of  Sp  with  respect  to  S* 
is  specified  by  an  (arbitrary)  analytical  ex- 
pression of  the  form  Sp  = A + B(S*)C  where  A , 

B , and  C are  constants. 

(b)  For  the  shelf  bathymetry  along  a particular  centrally 

located  isoline  of  £ , the  values  of  n are  de- 

termined which  will  yield  a constant  change  in  the 
time,  AT*  , required  for  a long  wave  traveling  at 
the  local  free  wave  celerity  to  proceed  from  seaward 
boundary  to  the  coast.  The  long  wave  travel  time, 

T , is  calculated  by: 


T = 

S 

n 

where  Sn  is  the  distance  along  the  isoline  of  £ , D0  is  the 
local  depth  relative  to  mean  sea  level,  and  g is  the  acceleration 
due  to  gravity.  The  relation  between  n and  T*  is  given  by: 

T*  = T*(T(Sn(n)))  , (39) 

The  incremental  values  of  T*  are  determined  from  equation  (38) 
subject  to  the  (arbitrary)  expansion  relationship  of  T*(T).  This 
relationship  is  a convenience  which  permits  an  additional  degree  of 
freedom  in  adjusting  the  relative  spacing  between  isolines  of  n • 
However,  if  T*  = T , the  value  of  AT*  is  that  which  divides  the 
total  long  wave  travel  time  by  an  integer  number  of  lines  of  n . 

The  selection  of  the  final  AT*  is  based  upon  a compromise  for  pro- 
viding adequate  resolution  for  the  wind  field  and  the  resulting  surge 
with  a minimum  number  of  points.  The  form  chosen  for  the  expansion 
is  T = A + B(T*)C  where  A , B , and  C are  constants. 

The  stretched  shelf  coordinate  system  provides  a grid  system 
with  a finer  resolution  near  the  coast  than  at  the  180-meter  iso- 
bath. The  expansion  curve,  S*(Sp)  , stretches  the  horizontal  reach 
of  the  grid  while  maintaining  a finer  grid  in  the  area  of  landfall 
of  the  hurricane.  In  this  manner  an  economy  is  achieved  in  terms  of 
the  number  of  grid  points  required  by  the  surge  program.  However, 
because  the  preferred  expansion  curves  dictate  the  locations  (in 
x,y  space)  of  the  depths  required  for  the  surge  calculations,  the 
depth  field  must  be  redefined  for  different  pairs  of  stretching 
functions.  Application  of  this  procedure  has  been  accomplished  for 
the  simulation  of  the  storm  surges  caused  by  Hurricanes  Carla, 
Camille,  and  Gracie.  Figure  24  shows  the  transform-generated  coast- 
line and  seaward  boundary  curves  for  the  western  gulf  coast  area. 

Also  shown  is  the  shelf  coordinate  system  (£,n)  which  was  employed 
in  the  Hurricane  Carla  surge  simulation.  The  track  of  the  hurricane 
is  the  dashline.  The  grid  system  extends  over  approximately  750 
kilometers  of  coastline.  The  stretched  shelf  coordinate  system 
(Figure  25)  is  attained  through  the  functional  relationships  speci- 
fied in  Figure  26  for  the  S*  axis  and  Figure  27  for  the  T*  axis. 

The  shelf  coordinate  system,  the  stretched  shelf  coordinate  system, 
and  similar  transformation  relationships  are  shown  in  Figures  28 
through  31,  for  the  grid  used  to  simulate  the  storm  surge  caused  by 
Hurricane  Camille,  and  in  Figures  32  through  35  for  the  grid  asso- 
ciated with  the  Hurricane  Gracie  storm  surge. 
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gure  25.  Stretched  shelf  coordinate  system  (S*,T*)  for  Hurricane  Carla 
surge  simulation. 


Functional  relationships  transforming  t,  to  S*  for  Hurricane  Carla 
surge  simulation. 


Carla  surge  simulation 


Jilt  i i.  ~ j — y.  ■>  , . r , . . 

of  simulated  Mississippi  Delta  shown  by  heavy  line  (this  extern 
seaward  end  of  the  Mississippi  River  levee) . 
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Figure  29.  Stretched  shelf  coordinate  system  (S*,T*)  for  Hurricane 
Camille  surge  simulation. 
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Figure  31.  Functional  relationships  transforming  n to  T*  for  Hurricane 
Camille  surge  simulation. 


Functional  relationships  transforming  r to  S*  for  Hurricane 
Gracie  surge  simulation. 


ttWspt 


UOl 


2.  Storm  Surge  liquations. 


The  vertically  integrated  form  of  the  quasi -linear  long  wave 
equations  in  a Cartesian  system  are  well  known  (e.g.  , Welander, 
1961).  The  appropriate  forms  for  storm  surge  computation  in  the 
stretched  shelf  coordinate  system  are  as  follows: 


3QS*  c ^ gD  3 , 

at  f ^t  Fu  as*  ^ ” B TS*  " °S*  ’ 

8QX*  pi)  3 

at  + f ^S*  + Fv  3T*  *-H  " V = TT*  ' °T* 


and 


3H  + ^2  iz  !c?  (FQc*)  + ^ |f*  (FQT*)]  = o , 


at  f2  n 3S* 


(40) 

(41) 

(42) 


where  Q is  the  volume  transport  per  unit  width,  t is  the  wind 
stress  divided  by  the  water  density,  a is  the  bottom  resistance 
stress  divided  by  water  density,  f is  the  Coriolis  parameter,  D 
is  depth  of  water,  H is  the  sea  surface  elevation  relative  to  mean 
sea  level,  and  Hg  is  the  hydrostatic  elevation  of  the  sea  surface 
corresponding  to  the  atmospheric  pressure  anomaly.  The  term  F ’is 
a variable  scale  factor  associated  with  the  conformal  part  of  the 
orthogonal  curvilinear  coordinate  system.  Specifically, 


(43) 


where  x and  y are  given  by  the  transformation  equations.  The 
terms  p and  v , are  additional  scale  factors  that  relate  to  the 
transformation  of  £ to  S*  and  n to  T*  , respectively,  and 
are  given  by: 


and 


as 

as* 


as„ 

-?n_  H 2I_ 

as  3T  3T*  ’ 

n 


(44) 


(45) 


Furthermore,  F2nv  corresponds  to  the  Jacobian  of  the  ti'ansforma- 
tion  in  the  sense  that 


A = / / F2uv  dS*  dT*  , 
R* 


(46) 


V 
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where 


Q = [Qsf  + QTf]^  , f 54 ) 

and  K is  a nondimensional  drag  coefficient  taken  as  2.5  * 10'3. 
o 

5.  Numerical  Algorithm. 

The  numerical  analogs  of  equations  (40),  (41),  and  (42)  are  based 
upon  centered  difference  approximations  of  all  terms  (see  App.  F). 

The  algorithm  treats  the  time  dependency  explicitly  and  employs  a 
computing  lattice  as  shown  in  Figure  56  in  which  the  transports,  Q, 
and  Q-p*  are  computed  at  the  same  location  but  are  staggered  in 
time  and  space  with  respect  to  the  water  level  anomaly.  This  scheme 
facilitates  a simple  representation  of  the  Coriolis  and  bottom  stress 
terms  in  the  difference  equations.  The  surge  model  allows  for  vari- 
able bathymetry  and  arbitrary  coastlines  which  are  represented  as  a 
series  of  connecting  straight  line  segments  situated  along  lines  of 
constant  S*  or  T*  . For  the  simulation  of  the  Hurricane  Camille 
storm  surge,  the  arbitrary  coastline  feature  of  the  algorithm  is 
mandatory  to  delineate  the  delta.  However,  in  the  more  usual  appli- 
cation, the  coast  is  a straight  line  in  the  computing  grid. 

The  difference  equations  for  Qg*,  Q-p*,  and  H at  interior 
points  of  the  computing  grid  are  given  by: 


Qs*( 

i , j ,n+l ) = (Gj  G2  + f At  Gj)/ (G2  + 

(f  At)2)  , 

(55) 

Qt*( 

i,j,n+l)  = (G j G3  - f At  G2)/(Gp  + 

(f  At)2)  , 

(56) 

and 


H(i, j ,n+l)  = H(i, j ,n-l) 


At 

F(i,j)2 


F(i+l,j)  Qs*(i+l,j,n)  - F(i-1, j)  Qg*(i-l,j,n) 
p(i)  AS* 


F(i > j +1)  QT*(iJ  + l,n)  - F(i,j-1)  QT*(i,j-l,n) 
v'(j)  AT* 


(57) 


where  i , j , and  n indices  express  the  S*  , T*  , and  time 
coordinates,  respectively,  and  At  is  the  numerical  time  step.  The 
quantities  Gj  % G2  , and  G3  are  given  by: 


Gj  = 1 ♦ 2 Kq  At  Q(i, j ,n-l)/D2  , 


(58) 
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G2  = Qs* ( i » j * n- 1 ) + f At  QT*(i,j,n-l) 


- H ( i - 1 


. j ,n  ) + llj{(  i - 1 , j , n ) | + 2 At  i ^ , 


g3  * Qt* C i . j , n - 1 ) - f A Qs*(i,j,n-1) 


WjAfQ')  AT*  H(i*j+1’n)  ' "B(i’j  + 1’n) 


+ Hpfi.j-l.n)  + 2 At  t , 


l)  = h 


D(i+l,j,n)  + D(i-l,j,n)  ♦ P(i,j+l,n)  + D(i,j-l,n) 


In  order  to  maintain  numerical  stability,  the  time  step  must 
satisfy  the  condition: 


Fjjv  AS*  AT* 

l'gp,_  J 


= I (a  AS*)2  + {v  AT*)2]'5  , (63) 

and  the  minimum  value  of  the  right-hand  side  of  relation  (26)  is 
implied.  Thus,  a search  of  the  grid  is  required  for  a proper  selection 
of  At  . However,  the  conditions  along  the  seaward  boundary  usually 
limit  the  time  step  because  of  the  deeper  water. 

The  values  of  a , v , F , and  0 given  by  equations  (44),  (45), 
(46)  and  (49),  are  required  at  the  appropriate  computational  grid 
points.  A numerical  spline  under  tension  is  employed  to  interpolate 
C given  the  1M  values  of  S (S*)  . The  scale  factor,  u , is 
determined  by:  ^ 

-m  - — 'Is--—  «“•> 


i = 1.2 


*4 


where  Uo)  and  £(IM+1)  are  ascertained  by  expanding  the  computing 
grid  one  increment  of  AS*  to  the  left  and  right,  respectively.  In 
a similar  fashion,  the  .INI  values  of  r,  may  be  determined  anil  the 
scale  factor  v given  by: 


v(j) 


n ( j 1 ) - n(j-l) 
2 AT* 


(65) 


j 


2,3 


JM-  1 . 


Since  the  computing  grid  may  not  be  extended  landward  or  toward  deeper 
water,  the  approximation  at  i = 1 and  i = JM  is  made  that: 

v ( 1 ) = [n(2)  - n(l)]/AT*  , (66) 

and 

v f JM)  = [n(JM)  - n(JM-l)]/AT*  . (67) 

The  area  in  the  x,y  plane  of  each  grid  block,  including  those 
from  expanding  the  grid  in  the  £ direction,  is  determined  by  numer- 
ical procedure  from  the  following  relation: 


A(i.j)  = 


Vl  Ci 


n.  £ 

J i-1 


(1*)? 

us 

(68) 

1 ,2  • ■ 

• 1 M+ 1 , 

1,2  • • 

• .JM  . 

The  scale  factor  associated  with  the  conformal  mapping  is  determined 
from  equation  (46)  at  intersectional  points,  excluding  the  coast  and 
seaward  boundary,  by  the  expression: 

r- , ; : , _ a ( i , j ) + A ( i , j + 1 ) + A ( i + 1 , j ) + A ( i + 1 , j + 1 ) * 

M1,J)  ' 4 AS*  AT*  Uf  i ) v(j) • (69) 

i = l ,2  • • • IM  , 

j = 2,3  •••  JM- 1 , 

and  for  grid  points  along  j = 1 or  JM  by  the  approximation: 
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■ -x  A ( i , i ) ♦ A ( i * 1 j ) 

1’-'  " 2 AS*  AT*  i.'(i’)  v(jl  ’ 

i = 1 ,2  • • • JM  . 

The  angle , 6 , relating  the  orientation  of  the  stretched  shelf 
coordinate  system  to  the  x,y  plane  is  determined  at  grid  points, 
excluding  j = 1 and  JM  , by  the  following  smoothing  procedure: 

0(i>j)  = l 0 ( i - 1 . j ) + 9‘(i,j-l)  + 9'(i  + l,j)  + e+(i,j  + l)  ♦ 4RT(i,j)  , 

(71) 

i = 1,2  •••  IM  , 
j = 2,5  •••  JM- 1 , 

4- 

where  9'  is  given  by  equation  (44).  Along  the  sea  and  coastline 
boundaries  we  have 

9 ( i . 1 ) * l[9+(i-l ,1)  + 29+(i,2)  ♦ e'  (i  + 1 ,1)  ♦ 40J  ( i , 1 ) I , (72) 


1,2  IM  , 


0(i,JM)  = a e'(i-l.JM)  + 20  (i.JM-1)  + 9 (i+1 , JM)  + 4e‘(i,JM)  , 


i = 1 ,2  • ■ • IM 


4.  Boundary  Conditions. 


A wall  condition  is  employed  in  the  surge  simulation  along  the 
coast  while  the  surface  elevation  anomaly  is  placed  in  equilibrium 
with  the  atmospheric  pressure  along  the  seaward  boundary.  Vanishing 
normal  derivatives  of  transport  are  specified  on  the  lateral  open 
boundaries.  This  condition  is  used  by  Jel esnianski  (1965,  1966)  and 
Forristall  (1974),  although  it  may  not  be  the  most  desirable  (Reid, 
1975) . 

In  the  x,v  plane  the  coast  is  curved  making  the  wall  (coast) 
and  lateral  boundary  conditions  difficult  to  apply  in  a rectilinear 
grid  system.  However,  the  stretched  shelf  coordinate  system  repre- 
sents the  coast  as: 

Qt*  = 0 , (74) 

or,  the  analog  is  simply, 


/>*  '*  • 


Qt*  ( i ,.JM,n+l ) = 0 , 


(75) 


i = 1,3  • • • IM  . 

Thus,  total  reflection  is  assured  at  the  boundary.  The  flux,  Qg*, 
along  the  coastline  is  calculated  from  equation  (55)  with  D taken 
as  the  two  point  average  of  the  local  fluid  depths  along  the  coast. 
The  local  depths  ranged  from  1.0  to  2.3  meters,  depending  on  the  sea- 
ward bottom  slope  and  surrounding  elevations.  The  water  elevations 
along  the  coast  are  computed  from  the  continuity  relation  as  given 
by  equation  (57)  with  the  following  substitution: 

F ( i ,.JM+ 1 ) Qt*  ( i , JM+1  ,n  - 1 ) = -F(i.JM-l)  (i  ,.JM-1  ,n-l ) , (76) 


i = 2,4  • • • IM- 2 . 

This  is  an  artifact  consistent  with  total  reflection. 

For  the  simulation  of  the  Hurricane  Camille  storm  surge,  the 
normal  routine  of  the  surge  program  was  interrupted  at  those  grid 
points  representing  the  protruding  Mississippi  Delta.  Along  this 
part  of  the  coast  (the  solid  heavy  line  in  Figures  28  and  29),  the 
normal  flux  was  set  to  zero  and  the  tangential  flux  was  determined 
from  either  equation  (55)  or  equation  (56)  with  the  Coriolis  term 
vanishing.  The  continuity  relation  is  altered  depending  on  the 
orientation  of  the  boundary  to  be  consistent  with  total  reflection. 

The  open  deep  sea  boundary  condition  is: 

H(i , 1 ,n+l ) = HR(i , 1 ,n+l)  , (77) 

i = 2,4  • • • IM- 1 , 


and 


H(i,2,n+1)  = [h ( i - 1 , 1 ,n+l ) + H(i+l,l,n+l) 

+ H(i-1 , 3,n+l)  + H(i+1 ,3,n+l) ]/ 4 , (78) 

i = 3,5  • • • IM-2  . 

Specifying  the  seaward  boundary  in  this  sawtooth  form  obviates  the 
calculations  of  Qg*  and  Qj*  along  the  boundary  since  they  are  not 
required  for  computations  at  interior  points. 

The  lateral  open  boundary  condition  requires  the  normal  gradient 
of  the  S*-directed  transport  to  vanish. 
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3QS. 

9S^  = 0 • <79> 

This  condition  implies  that  along  the  left  side  of  the  grid, 

Qs*d.j,n+1)  = Qs*(3,j,n+1)  , (80) 

j = 3,5  • • • JM 

and  along  the  right  side, 

Qs„(IM,j,n*l)  = Qs*(IM-2,j,n+l)  , (81) 

j = 3,5  •••  JM  . 

The  S*-directcd  transports  for  even  j and  i=2  or  IM-1  are  deter- 
mined by  the  average  of  the  two  neighboring  interior  values.  Addi- 
tionally, the  T*-directed  transports  are  required  along  i=2  or  IM-1 
and  are  calculated  from  equation  (56)  where  D is  again  the  average 
of  the  local  fluid  depth  along  the  boundary. 

Special  computations  are  required  at  seaward  corner  points  (2,2) 
and  (3,1).  Simultaneous  equations  are  solved  for  Qj*(2,2,n+1)  and 
Qs»(3,l,n+1)  with  the  approximation  that  Q-j*  (3, 1 ,n  + l ) is  the  aver- 
age of  the  two  neighboring  interior  values.  Similar  expressions  are 
used  at  the  right-hand  seaward  corner  points. 

Other  conditions  at  the  seaward  and  lateral  boundaries  were 
experimented  with  for  the  Hurricane  Carla  simulations.  These  in- 
cluded radiational  conditions  of  the  type  discussed  by  Reid  and 
Bodine  (1968)  and  Reid  (1975).  The  main  differences,  as  anticipated, 
were  close  to  the  lateral  boundaries  where  the  "flow  through"  condi- 
tion in  equation  (79)  gave  a more  realistic  rendition  of  the  long- 
shore flow. 

The  basic  numerical  model  was  used  successfully  by  Alvarez 
(1973).  Further  testing  of  the  algorithm,  where  in  some  cases  an 
analytical  solution  was  possible  for  comparison  purposes,  is  pre- 
sented in  Appendix  ('.. 

5.  Wj  njJ  and  Pressure  Fields. 

Of  major  importance  in  simulating  the  storm  surge  is  the  accurate 
portrayal  of  the  hurricane  wind  field  on  the  computing  grid.  More- 
over, the  representation  of  the  wind  must  be  time-dependent. 

The  Hydrometeorological  Section  of  the  National  Weather  Service 
(NWS)  provided  charts  of  the  surface  winds  at  10  meters  above  mean 
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sea  level  and  of  the  barometric  pressure  for  Hurricanes  Carla  and 
Camille.  One  method  for  representing  the  wind  field  is  to  digitize 
the  above  charts  for  the  hurricane  in  question.  Usually  the  wind  and 
pressure  data  depicting  the  storm  are  sampled  in  time  at  particular 
grid  points  and  interpolated  to  provide  the  necessary  input  to  the 
surge  model.  However,  such  input  data  are  tedious  and  laborious  to 
obtain  and  do  not  guarantee  better  quality  of  the  input  than  can  be 
obtained  by  analytical  representations  of  the  forcing  fields.  It  is 
this  latter  approach  which  is  used  in  this  study. 

For  a given  hurricane,  the  following  parameters  are  sampled  in 
time  from  the  NWS  charts:  position,  forward  speed,  central  pressure, 

radius  vector  (relative  to  the  storm  movement)  to  maximum  winds,  and 
the  maximum  winds.  With  an  analytical  representation  of  the  surface 
wind  and  pressure  fields  the  above  parameters  need  only  to  be  inter- 
polated in  time. 


The  wind  field  representation  prior  to  consideration  of  the  land 
influence  as  given  by  Jelesnianski  (1965)  is  employed  in  this  study. 
The  x,y-wind  components  for  a stationary  storm  are: 


x r, 


[ - (x-xh)  sin  <t>  - ( y-yh ) cos  ] 


and 


VR 

V = — 

y rh 


p(rh) 


[ + (x-x^  c°s  $ - (y-yh)  sin  ^ ^ 


where 


rh  = t(x‘xh)?  + (y_yh)2  ^ * 


(82) 


’ (83) 


VD  is  the  maximum  wind,  d>  is  the  ingress  angle,  and  F(r,  ) is 
3/2  1/2  h 

(rh/Rh)  if  rh  < Pn  or  ^Rh/rV  if  rh  — Rh'  Distance  from 

the  storm  center  (x^,  y ) to  the  region  of  maximum  winds  is 

The  translation  of  the  storm  p.  ''ides  an  alteration  in  the  wind  field 
which  is  carried  out  following  the  method  of  Jelesnianski  (1965). 

This  involves  the  vector  addition  to  the  above  field  of  a supplemen- 
tal velocity  whose  direction  is  parallel  to  that  of  the  storm  and 
whose  magnitude  depends  only  upon  r^  . The  resulting  modified 
values  of  and  V are  used  to  compute  the  wind-stress 

componeTts:  y 
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and 


t = K (V  2 ♦ V 2)  V , 
x x y x 


t = K (V  2 + V 2)  V . 

y x y y 


(84) 


This  constitutes  the  symmetric  wind-stress  field.  The  stress  compo- 
nents in  the  stretched  shelf  coordinate  system  are  determined  by 
applying  equation  (84)  in  equations  (47)  and  (48). 


The  surface  atmospheric  pressure  field  associated  with  the  hurri- 
cane is  given  by: 


P = P + (P 
o 00 


V 


(85) 


where  PQ  is  the  central  pressure  and  P^  is  the  far  field  pressure. 
The  term  Hg  in  equations  (40)  and  (41)  is  equivalent  to 

(P.  - p)/Pag- 

The  symmetric  surface  stress  field  does  not  reflect  the  influ- 
ence of  land.  Without  modeling  this  influence  the  analytically 
determined  wind  field  near  the  coast  would  not  be  consistent  with  that 
from  the  NWS.  A systematic  procedure  (App.  H)  is  employed  to  alter 
the  symmetric  wind  field  such  that  it  portrays  the  effect  of  land. 
These  winds  are  referred  to  as  deformed.  Mr.  Thomas  Reid,  Texas  AfjM 
University,  Department  of  Oceanography,  is  responsible  for  the  basic 
development  of  this  wind  model  (unpublished  manuscript).  The  pressure 
field  was  not  altered  in  the  nearshore  region  to  conform  with  the 
wind  field. 


Figure  37(a)  represents  the  symmetric  Hurricane  Carla  wind  field 
at  1600  Greenwich  mean  time  (G.m.t.),  11  September  1961  (approxi- 
mately 4 hours  before  the  storm  crossed  the  coast).  Isovels  are 
shown  in  meters  per  second.  The  analytically  deformed  wind  field  at 
this  time  is  shown  in  Figure  37(b).  The  influence  of  land  on  the  wind 
field  is  evidence  even  at  a time  when  the  storm  is  approximately  55 
kilometers  from  the  coast.  Figures  38  and  39  show  the  symmetric  (a) 
and  deformed  (b)  winds  at  2000  G.m.t., 11  September,  and  0000  G.m.t., 

12  September,  respectively.  The  rotation  and  reduction  of  the  wind 
vector  in  the  nearshore  region  to  reflect  the  land  influence  is 
illustrated  in  comparing  Figure  40(a)  for  the  symmetric  case  with 
Figure  40(b)  for  the  deformed  wind.  These  figures  correspond  to  the 
contours  shown  in  Figures  38(a)  and  38(b),  respectively.  The  arrows 
represent  the  wind  vector  placed  such  that  the  tail  is  at  the  com- 
putational point.  The  pressure  fields,  as  determined  from  equation 
(85) , for  Hurricane  Carla  at  the  above  times  are  shown  in  Figure 
41(a,b,c)  where  the  contours  are  in  millibars.  Although  the  pattern 
of  isovels  reflecting  the  land  influence  is  in  good  agreement  with 
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Figure  40.  Hurricane  Carla  winds  at  selected  points  at 
2000  G.m.t.,  11  September  1961  for  the 
symmetric  (a)  and  deformed  (b)  representations. 
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Figure  41.  Hurricane  Carla  atmospheric  pressure  on  11  Sep- 
tember 1961  at  1600  (a),  at  2"00  lb),  and  0000 
G.m.t.,  12  September  (c) ; contours  in  millibars. 
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those  provided  by  NWS,  the  problem  of  specifying  the  proper  deforma- 
tion is  not  completely  resolved. 

The  observed  high  winds  which  remain  along  the  coast  after  the 
storm  has  proceeded  inland  are  attained  by  setting  Rh  to  be  the 
distance  the  storm  center  is  from  that  coastal  area  of  persistent 
high  winds.  This  procedure  was  also  followed  in  modeling  the  winds 
from  Hurricane  Camille.  However,  in  the  pressure  expression 

was  not  allowed  to  increase  in  this  manner. 

The  hydrographs  at  selected  grid  points  along  the  coast  are 
presented  in  Figure  42  illustrating  the  improved  agreement  with 
respect  to  the  observed  water  levels  in  using  a deformed  wind 
forcing  function  rather  than  a symmetric  one.  The  observed  water 
levels,  corrected  for  the  astronomical  tide,  are  indicated  by 
squares.  The  computed  water  levels  have  been  raised  for  comparison 
purposes  to  correspond  to  the  local  sea  level  at  the  start  of  the 
computations  (1800  G.m.t.,  9 September).  The  solid  line  represents 
the  computed  surge  employing  the  deformed  wind  and  the  dashline 
shows  the  results  with  a symmetric  wind.  The  surge  computations 
were  formed  using  identical  boundary  conditions. 

The  symmetric  (a)  and  deformed  (b)  winds  for  Hurricane  Camille  at 
0000  G.m.t.  18  August  1969  (approximately  4 hours  before  the  storm 
crossed  the  coast),  0400  G.m.t.  and  0800  G.m.t.  on  the  same  date 
are  shown  in  Figures  43  through  45.  The  rotation  and  reduction  of 
the  wind  stress  for  this  storm  are  different  from  that  of  Hurricane- 
Carla  because  of  the  differences  in  the  coastal  configuration. 

Figure  46  fa,b)  shows  the  wind  vectors  for  the  symmetric  and  deformed 
representation  of  Hurricane  Camille  at  the  time  the  storm  proceeded 
inland.  The  pressure  fields  at  the  above  times  are  presented  in 
Figure  47|a,b,c).  The  wind  fields  associated  with  Hurricane  Gracie 
were  not  deformed  because  the  documentation  was  not  amenable  to  the 
deformation  procedure.  Also,  R^  was  not  increased  in  the  manner 
specified  above  after  the  storm  crossed  the  coast.  The  symmetric 
winds  at  1200  G.m.t.,  29  September  1959,  1600  G.m.t.  (the  time  the 
storm  proceeded  inland),  and  2000  G.m.t.  on  the  same  date  are  shown 
in  Figure  48.  Figure  49  presents  the  pressure  fields  at  the  above 
times. 
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figure  46.  Hurricane  Camille  winds  at  selected  points  at  0400 
G.m.t.,  18  August  1969  for  the  symmetric  fa)  and 
deformed  (h ) representations. 


Time  = 24  0 (hrs  ) 


Figure  48.  Hurricane  Gracie  symmetric  winds  on  29  September 
1959  at  1200  (a),  at  1600  (b) , and  at  2000  G.m.t. 
(c) ; isovels  in  meters  per  second. 
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Time  = 24.0  (hrs.) 


Time.=  28.0  (hrs.) 


§ Time  = 32.0  (hrs.) 
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Figure  49.  Hurricane  Gracie  atmospheric  pressure  on  29  September 
1959  at  1200  (a),  at  1600  (b),  and  at  2000  G.m.t. 

(c) ; contours  in  millibars. 


IV.  SIMULATION  OF  THE  FREE  WAVE  IN  AN  ANNULUS 


1 . Problem  Statement. 

If  the  wind  stress,  bottom  friction,  atmospheric  pressure,  and 
the  rotational  effect  of  the  earth  are  neglected,  then  the  equations 
in  Section  III  admit  a simple  analytic  solution  for  free  gravity 
waves  of  long  wavelength.  A study  using  curved  boundaries  will 
demonstrate  the  superiority  of  modeling  the  long  gravity  wave  in 
orthogonal  curvilinear  coordinates  over  those  models  which  employ 
rectilinear  coordinates.  Basically,  the  study  involves  comparing 
the  numerical  solution  of  a free-standing  long  wave  in  a 90°  section 
of  an  annulus  with  an  inner  radius  rj  of  393  kilometers  and  an 
outer  radius  T2  equal  to  2r^  (Figure  50).  The  annulus  is  bounded 
on  all  sides  by  a vertical  wall.  The  depth  of  the  basin  relative 
to  the  mean  water  level  is  assumed  everywhere  to  be  40  meters. 

The  analytical  solution  for  the  free-wave  oscillations  in  a 
section  of  an  annulus  may  be  obtained  with  some  modifications  from 
Lamb  (art.  191,  1932).  The  boundary  conditions  require  that: 

an 

— = 0 , at  r = ri  and  r = r2  , (86) 

3r 

and 

= 0 , at  0 = 0°  and  0 = n/2  . (87) 


The  analytical  solution  for  the  given  initial  conditions  is: 

°°  °°  r J ' (k  ri . 
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k is  the  mth  root  of 
m,n 


(88) 


(89) 


90 


\ 


f 


(90) 


Jn(krz)  Yn(kri)  ' Jn(krl)  Vkr2>  = 0 ' 

The  terms  .1  and  Y are  Bessel  functions  of  the  first  and  second 
kind,  respectively.  The  order  is  given  hy  tin  inter1',  subscript  and 
the  superscript  refers  to  differentiation  with  r.  spect  to  the  argu- 
ment. The  term  A in  equation  (88)  is  an  arbitrary  constant  repre- 
senting the  initial  amplitude.  The  first  a muthal  nodi  is  f r 
n = 2 and  the  lowest  value  of  k that  satisfies  equation  (90)  is 
approximately  1.340/rj  meters'1.  The  period  of  oscillation  for  this 
mode  (m  = 0 , n = 2)  is  25.85  hours. 

The  numerical  solution  of  the  free  wave  in  the  quarter  annulus 
is  sought  by  performing  the  integration  in  two  different  computing 
grids.  In  one  case  the  grid  is  rectilinear  (Figure  51),  and  in  the 
other  a polar  system  is  used.  In  the  rectilinear  system  the  outer 
and  inner  radii  of  the  annulus  (the  light  line  in  Figure  51)  is  simu- 
lated in  a stairstep  fashion.  Considering  a limitation  on  computer 
time  and  storage,  an  acceptable  rendition  of  the  curved  boundaries 
is  present  in  the  Cartesian  grid.  Proper  representation  of  the 
quarter-annulus  requires  a transport  point  at  corners  on  the  outer 
and  inner  boundaries.  Consequently,  the  rectilinear  boundary  is  not 
symmetric  about  tt/4.  The  locations  of  nine  hydrograph  positions  are 
indicated  by  small  boxes.  Although  the  computing  grid  is  43  by  43, 
only  1,052  points  are  used  to  represent  the  annulus.  In  this  grid 
system,  AS*  and  AT*  are  just  Ax  and  Ay  , respectively,  with 
Ax  = Ay  = 19.65  kilometers.  The  maximum  allowable  time  step  as 
determined  from  equation  (62)  is  approximately  700  seconds  where  F = 
u = v = 1.  Finally,  the  analog  form  of  the  long  wave  equations  is 
obtained  by  setting  those  terms  in  equations  (55)  through  (57)  to 
zero  which  are  neglected,  and  setting  F , p , and  v to  one.  The 
transports,  Q<.*  , Qy*  , are  now  aliases  for  Q and  Q , 
respectively.  x ^ 

The  numerical  algorithm  for  the  rectilinear  grid  system  must  be 
capable  of  determining  if  a computational  point  is  interior, 
exterior,  or  on  the  boundary  of  the  annu’us.  For  exterior  points,  no 
calculations  are  performed.  Grid  points  on  the  stairstep  boundary 
require  special  attention  subject  to  the  condition  of  a wall. 
Furthermore,  the  numerical  program  must  identify  and  apply  the  appro- 
priate wall  condition  depending  on  whether  the  point  in  question  is 
located  at  a boundary  comer  (an  H computation)  or  on  a segment  of 
the  boundary.  Clearly,  extensive  programing  and  computer  time  is 
required  to  accomplish  this  task. 

The  other  computing  mesh  which  is  used  for  obtaining  the  numerical 
solution  of  the  free  wave  in  the  annulus  is  the  polar  (or  stretched 
shelf  coordinate)  grid  system  (Figure  52).  The  transformation 
to  the  computing  grid  as  shown  in  Figure  53  is  accomplished 
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Figure  51.  Rectilinear  grid  representing  the  annulus 
location  of  hydrographs  shown  by  (□). 
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i - 1,2  •••  45  , 

j - 2,3  •••  22  , 

where  AS*  = 21.05  kilometers  and  AT*  = 19.65  kilometers.  Approxi- 

mately the  same  number  of  computational  points  are  employed  in  this 
system  as  that  of  Figure  51.  The  actual  area  of  the  annulus  is  well 
represented  in  both  computing  grids.  Although  the  inner  and  outer 
boundaries  are  straight-line  segments  in  the  numerical  grid  sense, 
the  portrayal  of  the  annulus  is  certainly  more  representative  than 
that  shown  in  Figure  51.  The  positions  of  nine  hydrographs  in 
Figure  52  (denoted  by  large  dots)  correspond  to  the  hydrograph  posi- 
tions in  the  rectilinear  grid  system.  For  the  polar  system,  the 
maximum  numerical  time  step  is  approximately  600  seconds.  As  a 
practical  note,  the  stretched  shelf  coordinate  system,  although  re- 
quiring less  sophisticated  programing  and  allowing  better  physical 
and  mathematical  portrayal  ot  the  boundaries,  generally  results  in 
a smaller  numerical  time  step  than  a rectilinear  system  of  comparable 
scale. 

2.  Results. 

For  comparison  purposes,  a time  step  of  540  seconds  is  used  in 
both  grid  systems  for  the  numerical  solution.  As  noted  above,  the 
polar  grid  program  should  compile  and  execute  faster  on  the  computer 
than  the  rectilinear  grid  program.  On  the  IBM  360/6S,  the  polar 
grid  program  required  0.47  minute  to  compile  and  0.39  minute  per 
100  iterations.  For  the  rectilinear  system,  the  program  consumed 
0.59  minute  to  compile  and  0.61  minute  per  100  iterations. 

The  initial  surface  elevation  condition  which  is  imposed  on  both 
grid  systems  is  that  3H/30  is  a constant.  For  the  rectilinear  grid 
system, 

H(i,j,0)  = .1  - ^tan"1  , (92) 


where  only  the  i,j  indices  appropriate  for  the  annulus  are  employed 
in  the  above  relation.  For  the  polar  grid  system,  the  same  initial 
surface  condition  is  determined  as: 
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H(i , j ,0)  - .1  ♦ .2  (i  - 1)  , (93) 

i = 1,2  •••  45  , 
j = 2,3  •••  22  . 

The  initial  condition  of  constant  azimuthal  slope  and  no  radial  vari- 
ation of  H along  a line  of  constant  0 implies  that  there  are 
several  modes  of  seiching  present;  however,  the  dominate  mode  is  for 
m = 0 and  n = 2.  This  is  indicated  by  Figures  54  and  55  which  show 
the  simulated  surface  topography  at  27  hours  in  the  rectilinear  and 
polar  grid  systems,  respectively.  The  solution  in  the  polar  grid 
is  more  representative  of  the  analytical  solution  than  that  shown 
in  the  Cartesian  system. 


The  hydrograph  at  9 = 0°  and  r = r;  is  shown  in  Figure  56  (a) 
where  the  water  elevation  as  determined  from  the  polar  grid  is  the 
solid  line  and  that  from  the  rectilinear  mesh  is  the  dashline. 

Figure  56  (b,c)  shows  the  hydrographs  along  6=0°  for  r = r 
(the  average  radius  of  the  annulus)  and  r = T2 , respectively.  The 
three  hydrographs  along  0 = tt/2  and  the  same  r positions  are  shown 
in  Figure  57  (a,b,c)  and  along  the  nodal  line  by  Figure  58  (a,b,c). 
The  average  period  of  oscillation  as  determined  from  the  hydrographs 
of  Figures  56  and  57  is  approximately  26  hours  in  the  polar  system 
and  28  hours  in  the  rectilinear  grid.  The  error  in  the  period  of 
oscillation  (about  8 percent)  for  the  rectilinear  system  is  most 
evident  in  the  figures  by  noting  the  lag  of  the  dashline  with 
respect  to  the  solid  line.  The  longer  period  of  oscillation  is 
directly  related  to  the  stairstep  boundary.  Effectively,  the  length 
of  the  basin  is  increased  by  the  reflections  introduced  by  these 
boundaries.  This  distortion  is  more  than  academic  since  many  recti- 
linear grid  models  of  enclosed  irregular  bays  require  adjustments 
to  reproduce  the  fundamental  seiching  mode. 

The  analytical  solution  at  any  point  is  a smooth  function  of 
time.  The  solid  lines  in  the  various  hydrographs  portray  this 
feature  better  than  the  dashlines  which  are  contaminated  by  high- 
frequency  spurious  oscillations.  The  nodal  line  in  the  polar  grid 
f . solution  as  evidenced  by  Figure  58  (a,b,c)  remains  fixed  at  0 = ir/4 

which  agrees  with  the  theory.  This  is  not  observed  in  the  recti- 
linear system.  Although  the  hydrographs  in  the  rectilinear  grid 
■ . are  not  positioned  exactly  on  0 = it/4  (actually,  about  44°),  the 

magnitude  of  the  oscillation  about  tt/4  is  approximately  2 to  4 
times  larger  than  expected. 
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Hydrographs  along  0 = r./2  and  r = rj 
(a),  rm  (b)  and  r2  (c)  for  the  polar 
grid  (solid)  and  rectilinear  grid 
(dashed) . 
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V.  STORM  SURGE  SIMULATION 


1.  Hurricane  Carla. 


Hurricane  Carla  was  an  immense,  slow-moving,  and  meandering  storm 
which  struck  the  gulf  coast  south  of  Port  Lavaca,  Texas,  at  2000  G.m.t. 
on  11  September  1961.  The  radius  to  maximum  winds  was  50  kilometers. 
The  atmospheric  pressure  drop  across  the  storm  was  75  millibars.  The 
maximum  sustained  winds  were  approximately  51  meters  per  second. 

For  numerical  stability  a time  step  of  180  seconds  was  used  in  the 
computations.  The  surge  simulation  was  performed  for  a 66-hour  period 
from  1800  G.m.t.  on  September  to  1200  G.m.t  on  12  September.  The 
boundary  conditions  are  those  given  in  Section  III.  The  wind  fields 
are  analytically  deformed  in  a manner  previously  discussed.  The  ob- 
served (squares)  and  computed  (solid  line)  water  levels  at  selected 
grid  points  along  the  coast  are  shown  in  Figure  59.  The  observed 
water  levels  have  been  corrected  for  the  astronomical  tide.  At  the 
start  of  the  computations,  the  water  surface  along  the  coast  from 
Sabine  Pass  to  the  Matagorda  Bay  area  was  elevated  by  approximately 
1 meter.  Consequently,  for  comparison  the  computed  water  levels 
have  been  raised  to  correspond  to  the  local  sea  level  condition.  The 
water  surface  topography  in  0.5-meter  contour  increments  at  46,  50, 
and  54  hours  are  shown  in  Figure  60.  The  computed  maximum  coastal 
surge  profile  corrected  for  the  astronomical  tide  is  presented  in 
Figure  61.  Indicated  on  the  figure  are  the  observed  high-crest  values 
from  tide  gages  and  water  marks.  The  depression  of  the  water  surface 
to  the  left  of  the  storm  track  (as  viewed  from  the  sea)  and  the 
low-peak  surge  in  the  Aransas  Pass  area  results  from  the  offshore 
wind  and  from  simulating  the  coast  as  a wall.  It  should  be  emphasized 
that  the  surge  model  does  not  consider  the  coastal  flooding  from  the 
shelf  surge  nor  includes  the  influence  of  the  increased  communication 
between  the  semienclosed  bays  and  the  open  sea.  While  the  computations 
agree  fairly  well  with  those  observed  for  the  peak  coastal  surges, 
especially  for  the  area  to  the  right  of  storm  landfall,  and  for  the 
Galveston  hydrography,  there  is  some  discrepancy  for  regions  far  from 
the  storm  track  such  as  Sabine  Pass.  This  might  be  attributed  to  the 
problem  of  proper  deformation  of  the  wind  field  (both  speed  and 
direction)  to  reflect  the  influence  of  land. 

Considerable  erosion  of  beaches  and  adjacent  offshore  areas  along 
the  Texas  coast  occurred  during  Hurricane  Carla  (U.S.  Army  Engineer 
District,  Galveston  1962).  Figure  62  shows  the  water  velocity 
V = Qs*i  * Qt* J ) /D  in  centimeters  per  second  at  selected  grid  points 
for  the  same  times  as  that  in  Figure  60.  Of  practical  importance  in 
assessing  the  erosion  potential  and  transport  of  the  material  is  the 
alongshore  current,  Qg*/ff  . In  shallow  water,  the  current  is  uniform 
from  the  surface  to  the  sea  bed.  Values  of  the  current  in  centimeters 
per  second  as  a function  of  time  at  selected  grid  points  are  presented 
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Figure  59.  Observed  (squares)  and  computed  (solid) 
water  levels  for  Hurricane  Carla. 
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Figure  60.  Computed  water  surface  topography  in  0. 5-meter 
contour  increments  for  Hurricane  Carla  at  46  (a} 
50  (b) , and  54  hours  (cl. 


astronomical  tide.  The  numbers  distributed  about  the  bottom  curve  representing 
the  coastline  are  observed  high  water  crests  in  meters  from  still  high  water  ma 

and  tide  gages  (underlined). 
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in  Figure  63.  Also  given  is  total  transport  in  106  cubic  meters  per 
second  as  a function  of  time  across  the  transect  lines  indicated  in 
the  figure. 

2 . Hurricane  Camille. 

Hurricane  Camille  traversed  the  extensive  low-lying  marsh  area 
of  the  Mississippi  Delta,  moving  over  the  shallow  waters  of  Breton 
and  Chandeleur  Sounds  for  4 hours  prior  to  landfall  near  Bay  St. 

Louis,  Mississippi,  at  approximately  0400  G.m.t.,  18  August  1969.  The 
radius  to  maximum  winds  was  about  37  kilometers.  The  atmospheric 
pressure , drop  across  the  storm  was  90  millibars.  The  maximum  sus- 
tained winds  were  approximately  60  meters  per  second  with  the  highest 
wind  gusts  estimated  at  90  meters  per  second. 

The  surge  simulation  extended  over  a 48-hour  period  beginning  at 
0000  G.m.t.  on  17  August.  A time  step  of  60  seconds  was  selected  for 
numerical  stability.  The  boundary  conditions  are  the  same  as  those 
used  in  the  Carla  simulation.  The  Camille  winds  are  analytically 
deformed  as  shown  in  Figures  43  through  45.  The  bottom  drag  coeffi- 
cient was  increased  to  5 * 10" 3 at  computational  grid  points  in  the 
area  bounded  by  6 <_  i £ 11  and  14  j JM  to  simulate  the  larger 
resistance  associated  with  the  low  flooded  marsh  lying  northeast  of 
the  simulated  Mississippi  Delta  (see  Figures  28  and  65). 

Figure  64  shows  the  observed  and  computed  water  levels  at  se- 
lected grid  points  along  the  coast.  At  Pascagoula,  Mississippi,  the 
water  level  increased  from  2.2  to  4.5  meters  in  a 2-hour  period 
indicating  a strong  convergence.  In  fact,  the  simulated  hydrography 
at  the  approximate  location  of  the  maximum  surge  (Bay  St.  Louis) 
showed  a similar  convergence  with  the  sea  level  rising  from  2.9  to 
7.0  meters  in  the  same  time  period  (figure  not  presented).  The  water 
surface  topography  in  1 -meter  contour  increments  at  selected  times 
from  24  to  36  hours  illustrating  the  surge  development  is  shown  in 
Figures  65  through  73.  The  effect  of  the  delta,  modeled  as  a wall 
protruding  from  the  coast,  is  visible  in  these  figures  by  the  presence 
of  high  water  along  the  eastern  side  of  the  wall.  The  track  of  the 
storm  with  respect  to  the  delta  is  such  that  the  circulation  is 
inhibited  in  this  area  resulting  in  the  setup.  The  water  velocity, 

V , at  the  selected  times  and  grid  points  as  shown  in  Figure  74  are 
in  accord  with  these  observations.  The  computed  maximum  coastal 
surge  envelope  east  of  the  delta  presented  in  Figure  75  shows  that 
the  calculated  peak  s, ”ges  are  approximately  that  observed  in  the 
Gulfport  to  Pascagoula  area  (U.S.  Army  Engineer  District,  New  Orleans, 
1970).  A test  run  employing  a constant  wind-stress  coefficient,  K , 
of  3 x 10"6  with  the  same  boundary  conditions  produced  a similar 
maximum  surge  envelope.  There  is  qualitative  agreement  between  the 
observed  and  predicted  surge  envelopes.  This  is  shown  by  the  rela- 
tively low  water  elevation  of  2.7  meters  near  Hopedale,  Louisiana, 
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Figure  63.  Values  of  alongshore  current,  Qs./D  , in  centimeters  per  second  as  a function  of 
time  at  selected  grid  points.  Total  alongshore  transport  in  cubic  meters  per 
second  is  indicated  by  I . 


Computed  water  surface  topography  in  1.0-meter  contour  increments 
in  the  area  of  Breton,  Chandeleur,  and  Mississippi  Sounds  for 
Hurricane  Camille  at  24  hours. 


Pascagoula 


Computed  water  surface  topography  in  1.0-meter  contour  increments  in 
the  area  of  Breton,  Chandeleur,  and  Mississippi  Sounds  for  Hurricane 
Camille  at  26  hours. 


Computed  water  surface  topography  in  1.0-meter  contour  increments 
in  the  area  of  Breton,  Chandeleur,  and  Mississippi  Sounds  for 
Hurricane  Camille  at  27  hours. 
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Figure  72.  Computed  water  surface  topography  in 
Breton,  Chandeleur,  and  Mississippi  : 


nomical  tide.  The  numbers  distributed  about  the  bottom 
curve  representing  the  coastline  are  observed  still  high 
water  crests  in  meters. 


in  the  marsh  area  south  of  Lake  Borgne,  and  a high  water  mark  of  4.8 
meters  along  the  eastern  side  of  the  Mississippi  River  levee.  A 
part  of  the  discrepancy  between  the  observed  and  computed  surges  re- 
sults from  not  considering  the  flooding  of  low  coastal  areas  and  the 
open  communication  of  Mississippi  Sound  with  Lake  Borgne  which  is, 
in  turn,  connected  with  Lake  Pontchartrain . The  tide  gage  at  Shell 
Beach  located  on  the  southern  end  of  Lake  Borgne  crested  at  3.4 
meters.  Another  test  simulation  was  conducted  with  the  hurricane 
track  moved  18  kilometers  to  the  east.  The  results  showed  greater 
seaward  circulation  in  Chandeleur  Sound  reducing  the  setup  along  the 
levee  by  approximately  1.5  meters.  The  location  of  the  peak  surge 
was  shifted  to  the  east  but  the  surge  magnitude  was  not  reduced. 

3.  Hurricane  Gracie. 

Hurricane  Gracie  crossed  the  coast  near  Beaufort,  South  Carolina, 
at  1600  G.m.t.  29  September  1959.  The  radius  to  maximum  winds  was 
25  kilometers.  The  atmospheric  pressure  in  the  eye  was  950  millibars. 
The  maximum  sustained  winds  were  approximately  50  meters  per  second. 

A Cartesian  grid  for  simulating  the  surge  caused  by  Gracie  is 
shown  in  Figure  76.  The  reach  and  grid  spacing  of  this  system  are 
comparable  to  the  curvilinear  grid  shown  in  Figure  32.  The  thin 
dashlines  are  a part  of  the  transform-generated  coast  and  seaward 
boundary  curve  (180-meter  depth  contour) . The  surge  computations  are 
performed  with  both  grid  systems  using  the  same  boundary  conditions 
and  forcing  functions.  The  stairstep  boundary  representing  the  coast 
is  the  heavy  solid  line  in  Figure  76.  Another  difference  between  the 
models  other  than  the  coastline  representation  is  that  the  Cartesian 
system's  boundary,  j = 1,  extends  farther  seaward  into  deeper  water 
than  that  of  the  curvilinear  system. 

For  numerical  stability  a time  step  of  150  seconds  was  used  in 
the  computations  with  the  curvilinear  grid.  The  time  step  in  the 
rectilinear  system  was  taken  as  75  seconds.  The  surge  simulation  was 
performed  for  a 48-hour  period  with  zero  hour  corresponding  to  1200 
G.m.t.  on  28  September. 

The  observed  and  computed  water  levels  along  the  coast  in  the 
curvilinear  (solid)  and  rectilinear  (dashed)  systems  are  shown  in 
Figure  77.  The  inserts  on  Figure  76  show  the  simulated  hydrograph 
locations.  Figures  78  and  79  present  the  water  surface  topography  in 
the  curvilinear  and  rectilinear  grids,  respectively,  at  the  same  times 
and  for  the  same  contour  increments.  The  depression  of  the  water 
surface  to  the  left  of  the  storm  track  exhibited  in  both  figures  is 
similar  to  that  from  the  Carla  simulation.  At  the  Savannah  River 
entrance,  the  observed  water  level  corrected  for  the  astronomical 
tide  is  depressed  after  27  hours  reflecting  the  offshore-directed 
wind.  The  effect  on  the  coastal  surges  and  water  surface  topographies 
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in  the  rectilinear  (□)  and  curvilinear  (•  ) grids.  Dots  on  the  dashlines  are  comer 
points  of  the  curvilinear  grid  given  in  Figure  32. 
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Figure  77.  Observed  (squares)  and  computed  water  levels 

for  Hurricane  Grade  in  the  curvilinear  (solid) 
and  rectilinear  (dashline)  grid  systems. 
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resulting  from  the  stairstep  portrayal  of  the  coast  boundary  are 
shown  through  the  comparison  of  Figures  78  and  79.  The  curvilinear 
model  shows  the  maximum  coastal  water  level  at  the  time  of  landfall 
of  approximately  3 meters  while  the  Cartesian  system  predicts  4 
meters.  Furthermore,  the  peak  surge  of  4.5  meters  in  the  Cartesian 
system  occurred  at  the  corner  point  (22,  19)  over  2 hours  after 
storm  landfall. 

The  water  velocity,  V , from  the  curvilinear  system  is  presented 
in  Figure  80.  The  figure  also  shows  the  development  of  a cyclonic 
current  in  the  region  from  Cape  Romain  to  Cape  Fear  resulting  from 
the  curved  coastline  and  the  extensive  Frying  Pan  Shoals  off  Cape 
Fear.  A similar  current  pattern  is  not  observed  in  the  rectilinear 
system. 

The  computed  high  water  envelopes  from  both  grids  are  shown  in 
Figure  81.  The  stairstep  rectilinear  coastline  is  projected  upon  the 
curvilinear  coastline  providing  a common  basis  for  the  comparison. 
Both  systems  predict  a maximum  surge  comparable  to  that  of  the  near- 
est observed  high  water  mark.  The  effect  of  the  rectilinear  coast- 
line on  the  surge  envelope  is  seen  in  the  figure  by  the  relative 
water  level  maximums  at  the  interior  corner  grid  points.  The  pri- 
mary contrast  of  the  coastal  envelope  of  computed  surges  in  the  two 
different  grid  systems  is  the  greater  variability  of  that  from  the 
rectangular  grid.  This  leads  to  spurious  results  as  can  be  seen 
particular  for  Savannah  River. 
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Figure  80.  Computed  *atcr  velocity  in  centimeters  per  second  at 
selected  grid  points  for  Hurricane  Gracic  at  24  la), 
28  (b),  and  52  hours  (c)  in  the  curvilinear  grid. 

See  Figure  5 for  identification  of  capes  referred  to 
in  the  text. 
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VI.  SUMMARY  AND  CONCLUSIONS 


The  development  of  the  conformal  mapping  problem  was  completed 
for  transforming  the  interior  region  bounded  by  two  arbitrarily 
shaped  curves  and  two  parallel  lateral  boundaries  into  a rectangle 
in  the  image  plane.  Specifically,  the  coastline  and  seaward  boundary 
curve  taken  as  the  180-meter  depth  contour  are  transformed  into  the 
image  plane  as  constant  values  of  n . The  bicurve  fitting  procedure 
was  successively  tested  on  five  regions  of  the  Continental  Shelf  in 
the  Gulf  of  Mexico  and  east  coast  of  the  United  States.  The  distinct 
advantages  of  the  orthogonal  curvilinear  coordinate  system  over 
rectilinear  representations  of  the  coast  are  the  accuracy  and  ease 
of  specification  of  the  boundary  conditions.  For  long  wave  calcula- 
tions, the  numerical  simulation  of  a free  long  wavelength  gravity 
wave  in  a quarter  annulus  demonstrated  the  superiority  of  modeling 
the  boundaries  in  curvilinear  coordinates. 

A two-dimensional  time-dependent  numerical  storm  surge  model 
employing  orthogonal  curvilinear  coordinates  is  presented.  Since 
the  transformation  is  conformal,  the  associated  modification  of  the 
equations  is  minimized.  The  final  coordinates  allow  for  the  greatest 
resolution  near  the  coast  over  the  area  of  principal  storm  surge 
development  and  modification.  A second  transformation  is  performed 
independently  on  each  of  the  curvilinear  coordinates  to  provide  the 
desired  spatial  resolution  in  the  shallower  region  near  hurricane 
1 and fa 1 1 . 

Results  of  the  numerical  simulation  in  transformed  coordinates 
of  the  storm  surges  along  the  gulf  coast  induced  by  Hurricanes  Carla, 
Camille,  and  Grade  which  crossed  the  east  coast  are  presented.  The 
surface  wind  fields  are  analytically  represented.  Additionally,  a 
systematic  procedure  is  employed  to  deform  the  symmetric  wind  fields 
of  Carla  and  Camille  in  the  nearshore  region.  The  deformed  Carla 
wind  representation  produces  better  agreement  between  the  computed 
surges  and  that  observed  than  a similar  simulation  with  a symmetric 
wind.  Although  the  pattern  of  windspeeds  produced  by  the  deformation 
procedure  is  in  good  agreement  with  those  provided  by  NWS,  the 
problem  of  specifying  the  proper  orientation  of  the  wind-stress 
vector  in  the  nearshore  region  requires  further  study. 

The  surge  simulations  indicate  that  the  model  produces  results 
in  good  agreement  with  the  observed  peak  surges  and  hydrographs, 
especially  in  the  area  to  the  right  of  landfall.  Moreover,  the 
results  from  the  simulation  of  the  Hurricane  Gracie  surge  indicate 
superior  rendition  in  the  curvilinear  grid  compared  with  that  of  the 
Cartesian  grid.  Certainly  one  contributor  to  any  discrepancy  between 
computed  and  observed  water  levels  may  result  from  the  portrayal  of 
the  coast  as  a wall.  Proper  inclusion  of  attached  bays,  lakes,  and 
flooding  of  low-lying  areas  should  be  pursued. 
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APPENDIX  A 


CONFORMAL  MAPPING  COEFFICIENTS 
FOR  THE 

GULF  COAST  REGION  FROM  ATCHAFALAYA  BAY  TO  APALACHEE  BAY 


B 

= 0.  1991  53  IE 

02 

B 

o 

= 0. 1 1 291 02E 

03 

X 

= 317.5 

n 

B 

C 

n 

n 

1 

-0. 1861 91 7 E 0 2 

— 0 .20 1 3675E  01 

2 

— 0. 2 8 04 65 2 E 02 

0.6239547E  01 

3 

0 • 64  03  1 3 1 E 01 

— 0.5958 1 48E  01 

A 

-0.2530149E  01 

0.5130100E  01 

5 

0.  153 1 1 2 IE  01 

0,  22 73 1 9 OE  Cl 

6 

0. 1498879 E-0 1 

0.351 2568E  01 

7 

-0.  49531  18E  00 

— 0.1 6 06539  E 01 

8 

— 0 • 1035305E  01 

0.1131 84  2E  01 

9 

-0 • 233444  6 E 00 

-0  .6252772E  — 01 

10 

0. 5907305E-01 

0.6906291 E 00 

1 1 

0.4296697E  00 

-0.1417943E  00 

12 

-0. 2997332E-0 1 

0.1 326714E-01 

13 

0. 4310032E- 0 1 

— C. 2609431 E 00 

1 4 

-0.6075356  E—  0 2 

0.1616244E— 01 

15 

0.111 7075E  00 

— 0.4899038E— 01 

16 

0 . 5209380E- 0 1 

0. 7586  825E- 0 1 

17 

0. 7516444E-02 

-0.2902573E-01 

18 

-0.  1 997  32  IE—  0 1 

0.  34  4 0 02 4E  — 01 

19 

— 0 . 579  306  2 E—0 3 

— 0.339441 5E—  0 1 

20 

0.  13441  83E-01 

0.2085005E— 01 

21 

0. 1558170E-0 1 

— 0, 5 30  9 36 2E  — 02 

22 

0. 3357759E— 0 2 

0. 1572620C-01 

Note:  all  quantities  are  in  x ,y  units. 


133 


***. 


1 


n 

B 

C 

n 

n 

23 

0.453654  3E- 02 

- 0.4705735E  — 02 

24 

—0 .6767012  E— 0 2 

0.4449  1B9E- 02 

25 

0. 5650939E-03 

0.154 1801 E— 02 

26 

— 0.3992960E— 02 

0*691585 2E— 02 

27 

0. 2905252E— 02 

0.2670427E-02 

2B 

-0. 5924194E- 03 

0.  31  55482E-02 

29 

0 . 396  4 6 1 0 E— 0 3 

— 0.1662592E— 02 

30 

-0.  87095356-03 

0. 5 55  0 05  2 E — 03 

31 

0 .4148923E-0  3 

-0.1 224325E- C2 

32 

0. 4497680E-0  3 

0.1261 886E-02 

33 

0. 555896 CE—  0 3 

-0.265304  OE-03 

34 

-0  • 40540 87 E—  0 3 

0.4554733E-03 

35 

-0. 252 979 8 E—  03 

— 0.5726672E  — 03 

36 

— 0.281 328  8E- 0 3 

0. 505077  8E—  04 

37 

0 . 3256833E— 0 3 

-0. 1338352E-03 

38 

0.  1 094  8 75E— 03 

0.  1 940636E-03 

39 

0. 1 343705E—  0 3 

— 0.1790418E— 03 

40 

-0.  £ 795955 E—  04 

— 0 ,2  37  754  3E  — 04 

41 

-0 • 8443834F- 05 

-0.166 1 804E- 03 

42 

-0.161 1809E-04 

0.9 14 19806-04 

43 

0. 4455329E—  04 

0. 5 84  22  C9E  — 05 

44 

-0. 212978  3E—  0 4 

0.6875678E- 04 

45 

-0.  139392  0E—  05 

— 0.3022245E— 04 

46 

-0. 3849607E- 04 

0.  5 1 77958E- 05 

47 

0.  142  74  32  E— 0 4 

-0 .1 272952E-04 

48 

-0.  389254  IE— C 5 

0.  1 8866626-04 

49 

0. 1366723E-04 

-0. 151 4819E- 05 

50 

-0. 9772998E— 05 

0.111 4 4 08  E— 04 

51 

-0. 376284 2E- 05 

-0.741  07C8E-05 

52 

-0. 761 0149  E— 0 5 

0.6834680E-05 

53 

0. 851 1489E- 06 

-0.3682432E— 05 

54 

-0 .8219286E-06 

0.50198816-05 

55 

0. 2588446E—  05 

— 0.25  87  597  E— 05 

56 

-0.  164  524 2E-  05 

0.  794  1880E-06 

57 

0.9023091 E-06 

— 0.2933107E-05 

58 

-0.  1668636E— 05 

-0. 71 80398E-07 

59 

0 .8869575E- 06 

-0.91  1 4 828E  — 06 

60 

-0. 4799319E-06 

0.1 0 67656  E — 05 

61 

0. C261487E- 06 

-0.  4 895755E-  C6 

62 

—0 . 523  7882  E—  06 

0.2304607E— 06 

63 

-0.  47596566-07 

-0.89250436-06 

64 

-0. 25988156-06 

0.  191  2503E-06 

65 

0. 32  *394 1 E—  06 

— 0.2104  3 35E  — 06 

66 

O. 293C804E- 07 

0.  373  0255E-C6 

67 

0.  102003BE— 06 

-0.141 70306-06 

68 

-0.  24221  00E-06 

0.591681 16-07 

134 


n 

B 

C 

n 

n 

69 

-0. 

588  016  08- 0 8 

-0. 1 276744E- C6 

70 

-0. 

462  6 68  4 E — 0 7 

0.1424694F-06 

71 

0. 

1 153690E-06 

0.  1 71  9058E-07 

72 

-0  . 

1 69  222  0 E— 0 7 

0.888053QE- 07 

73 

0. 

6679772E-09 

-0.6668830E-07 

74 

-0. 

4858599F- 07 

0. 169  8 34  PE— C7 

75 

0. 

2245389E— 0 7 

-0.2062551E-07 

76 

0. 

4829538E- 0e 

0. 5085036E  — 07 

77 

0 . 

22 10  1 85  F—  0 7 

-0.2895261E-08 

78 

-0. 

1 637371 E—  0 7 

0.1 01 3160F-07 

79 

0. 

5 77  758  3E- 0 9 

-0.  1960430E— 07 

80 

-0  . 

6 394433  E—  0 8 

0 .8530  744E-08 

81 

0. 

1182490F- 07 

— 0.2038007E  — 08 

82 

0. 

650 139  3 E—  0 9 

0.9465726E- Ce 

83 

0. 

26085  79E—  08 

-0.609281 5E-08 

84 

-0. 

4322956E- Oe 

0.  1 981  51  OE-  08 

85 

0. 

1 574734  E-08 

-0  .3446277E  — 08 

86 

-0. 

47261 93E-1 0 

0.414  6546F  — 08 

87 

0 . 

24 1 5840E- 08 

— 0. 909  2 067C- 09 

88 

-0. 

1 000588E-08 

0.1491  1 53  E — 08 

89 

0. 

6 74  27  5 3E- 09 

-0.  1408612E-08 

90 

-0. 

672  8380  E— 0 9 

0 .1087402F- 08 

91 

0. 

1012167E-08 

-0. 1 87931 3E-09 

92 

-0  . 

301 6696E- 0 9 

0.91631 99E—  09 

93 

0. 

323  686  0 E— 0 9 

— 0.3496718E  — 09 

94 

-0. 

3145593E- 09 

0.4997C83F— 09 

95 

0. 

335 3 65 3 E— 09 

— 0.1249788E—  09 

96 

-0. 

3262997F— 1 0 

0.440561 1 E— 09 

97 

0 . 

1631445E- 09 

— 0.1S01641E— 09 

98 

-0. 

13751 70E-09 

0. 1637698E-09 

99 

0. 

1 072999E- 09 

-0.  9262  971  E—  1 0 

1 00 

-0. 

1 384905F- 1 1 

0.1 895707E-09 

101 

0. 

1 08674  2 E—  0 9 

— 0.1663476E  — 10 

102 

-0. 

5080765E- 1 0 

0. 689067 OE— 10 

1 03 

0. 

1 376387E-1 0 

-0.751 0520E- 10 

104 

-0. 

2277  74  3E— 1 0 

0.4346501E— 1 0 

105 

0. 

5459 14  6 E— 1 0 

— 0.7598091E— 1 1 

106 

0. 

537051 7E- 1 1 

0.4664 877E— 10 

1 07 

0* 

1 35594  IE-  1 0 

- 0.2425636E- 10 

1 08 

-0. 

242  22 1 0 F—  1 0 

0. 4706750E— 1 1 

109 

0. 

7556705E- 1 1 

-0.  1434  774E-10 

1 10 

0. 

3368470  E— 1 2 

0.201 09S5E- 1 0 

1 1 1 

0. 

1 44  7254  E— 1 C 

0.1 050866  E — 1 1 

1 12 

-0. 

5194694E- 1 1 

0. 715031  IE- 1 1 

1 13 

-0. 

4 364 368 E- I 3 

— 0.9073677E— 1 1 

114 

-0. 

475  20  8 3E— l 1 

0.  31 38752E-1 1 

135 


T 


n 

B 

C 

n 

n 

X 15 

0 .4S53522E- 1 1 

-0.6452B75E- 12 

1 X 6 

0.  7602  034  E—  1 2 

0.6237037E— 1 1 

X X 7 

0.  1944S47E- 1 l 

-0.111 44 12E- 1 1 

1X8 

-0  *2436083  E— 1 1 

0 .1234543E-1 1 

X 19 

0. 191 7147E-1 2 

-0.  1740079E-1  1 

120 

-0.4928470E-12 

0.1649892E- 1 1 

121 

O.  1484346E-1  1 

0.5327387E— 1 3 

X 22 

-0.  1348819E-  12 

0.1  120676E— 11 

123 

0.3641300E-1 2 

-0 .5990049E— 12 

124 

-0. 5487838E-1 2 

0.2460370E— 12 

125 

0 • 24 0 6 55  8 E-  1 2 

— 0,344  8 14  3E—  12 

126 

-0. 791 1 185E-1 3 

0.4439149E— 12 

127 

0. 3179581E- 1 2 

— 0.5831951E— 13 

128 

-0.861 194  3E-1 3 

0.2390251E— 12 

1 29 

0.  72  3 84  7 0E-1  3 

— 0.1 814728E  — 12 

1 30 

-0 . 1 36 1 322E- 1 2 

0.557  39  03E— 13 

1 31 

0. 681 789  2 E— 1 3 

— 0.97457  36E  — 1 3 

132 

-0.  1 4 3696  9E- 1 2 

0. 1 1 52394E-12 

133 

0 . 8566666  E— 1 3 

-0.1 163015E- 13 

1 34 

-0. 3499625E-1 3 

0.5682  837 E - 1 3 

135 

-0.78356 16E- 1 5 

-0.551  0710E-  13 

1 36 

— 0 • 4394  79  2E— 1 3 

0.1780  100E-1  3 

1 37 

0. 23494 19E- 1 2 

— 0.  1 2 72559E  — 1 3 

138 

0 • 43958 1 0 E— 1 4 

0 .423  3778E— 1 3 

1 39 

0.  1 9 7682  BE— 1 3 

-0.381 6127E-14 

140 

-0. 1 68765 1 E- 13 

0.7731344E- 14 

141 

-0. 1897316E-14 

-0. 1672038E- 13 

142 

— 0. £78145  IE— 1 4 

0.  9221  1 eOE-14 

143 

0,  1 174166E- 1 3 

0 • 2324292E- 14 

144 

0. 6883079E-1 5 

0.1 088460E— 13 

145 

0. 211 6535E- 14 

- 0. 4 276968E  — 14 

1 46 

-0 . 55740  0 5 E— 1 4 

0.4950588E— 15 

147 

0.  13  2 1 55 IE— 1 4 

— 0.3220331E  — 14 

148 

0 .5106463E- 1 S 

0.424  1 152E-  14 

149 

0.  3148336E-1  4 

0.251 4837E-15 

150 

-0. 61591 13E- 1 5 

0.1642e51E-14 

151 

0 .450  7668  E— 1 6 

— 0.1803004E— 14 

152 

-0. 1 C82406E- 1 4 

0.54?0  030E  — 1 5 

153 

0 .919321 9E— 1 5 

-0.2248083E- 15 

1 54 

0. 21521 37E-1 6 

0.1 072704E-14 

1 55 

0.  544 053 9E-  1 5 

-0.201  94S9E-15 

1 56 

— 0.4207780  E- 1 5 

0.27  8 7050F— 15 

1 57 

0.  124  0553E— 1 5 

— 0.3249225E— 15 

158 

— 0.1 8099  34 E- 15 

0.2  89  6084E— 15 

159 

0. 21 3424 1 E— 1 5 

—0.65  P6600E— 16 

160 

— 0.641496  6E-  16 

0.  253  3 1 24E  — 1 5 

1 36 


r 


n 

B 

C 

n 

n 

1 15 

0.45535226- 1 1 

-0.6452875E - 12 

116 

0.  7602  034  E— 1 2 

0.6237037E— 1 1 

1 17 

0.  1944S47E-  1 1 

- 0. 1 1 1 44  1 2E-  1 1 

1 1 a 

— 0 .2436083  E— 1 1 

0.1234543E— 11 

119 

0.  191  7147E-12 

-0.  17  4 00  7 9E  — 1 1 

1 20 

— 0.4920470E- 1 2 

0.164  9 892E-  1 1 

121 

0.  1404  346E-1  l 

0.5327387E— 1 3 

1 22 

-0.  1348919E-  1 2 

0.1 120676E-11 

123 

0. 3641 300 E-l 2 

-0 .599  0 04  9E—  12 

1 24 

-0. 6487838E-1 2 

0.2460370E— 12 

125 

0. 2406S58E- 1 2 

-0. 3448  143E- 12 

126 

-0. 791 1 185E-1 3 

0.4439149E— 12 

1 27 

0. 3179581E- 1 2 

- 0. 583 1 951 E— 1 3 

128 

-0.061  194  3E-1 3 

0 .2390251E—  12 

129 

0.  723  84  7 0E—  1 3 

-0.1  81  4728E-12 

1 30 

-0. 1 36 1 32  2E-  1 2 

0.557  39  03E— 13 

1 31 

0. 681 7892E-1 3 

— 0.97457  36E  — 1 3 

132 

-0.  1 4 36969E-  1 3 

0. 1 1 52394E-12 

133 

0,8566666  E— 1 3 

-0.1 1630 15E- 13 

1 34 

-0. 3499625E-1 3 

0.5682837E— 13 

135 

-0.78356 16E- 1 5 

-0.551  0710E-  13 

1 36 

— 0 . 4394  792E— 1 3 

0.1780  1 00  E— 1 3 

1 37 

0. 23494 19E- 1 3 

— 0.  1 2 72599E  — 1 3 

138 

0 • 43958 1 0 E— 1 4 

0.42  3 3 77PE— 1 3 

1 39 

0. 1 976828E-1 3 

-0.381 6127F-14 

1 40 

-0. 168765  IE-  1 3 

0.77  31 344E—  14 

141 

-0. 1897316E-1 4 

-0. 1672038E-13 

142 

-0. 678145  IE- 1 4 

0.0221 1 e0E-14 

143 

0.  11  74 166 E-l 3 

0.23242 92E— 14 

144 

0. 6883079E— 1 5 

0.1 088460E-13 

145 

0. 21 1 6535E- 1 4 

- 0.4276968E- 14 

1 46 

-0 . 55740  05  E—  1 4 

0.4950588E— 15 

147 

0.  132155  IE— 1 4 

-0.3220331E  — 14 

148 

0.5106463E- 15 

0.424  1 152E-  14 

149 

0. 31 48336E-1 4 

0.251 4837E— 15 

150 

-0.815911 3E- 1 £ 

0.  1642e51E- 14 

151 

0 . 450  7668  E— 1 6 

-0.180  3004E- 1 4 

152 

-0.  1 C82406E- 1 « 

0.5420  C3  OE  — 1 5 

1 53 

0 .9193219E— 1 5 

— 0.224  8 083E— 15 

1 54 

0. 21521 37E-1 6 

0.1  07  2 704  E— 1 4 

1 55 

0. S440539E-  1 £ 

-0.201  9489E-15 

1 56 

— 0. 4207  780  E— 1 S 

0.2787050F- 15 

157 

0.  1240553E-1 5 

— 0.3249225E  — 15 

156 

-0. 1809934E- 15 

0.2896084E- 15 

1 59 

0. 2134241 E-l 5 

— 0 . 65  86600E- 16 

160 

-0. 641 4966E-  1 6 

0. 253312  4E  — 1 5 

137 


n IW1.WWI  <.  .. 

* 


#4 


APPENDIX  B 


CONFORMAL  MAPPING  COEFFICIENTS 
FOR  THF 

EAST  COAST  REGION  FROM  CAPE  KENNEDY  TO  PAMLICO  SOUND 


3 = 0.991  7248E  01 

B = 0.7052965E  02 
o 

X = 250. 0 


n 

B 

c 

n 

n 

1 

-0. 4282722E 

02 

-0.  62954  06E 

02 

2 

—0 « 282 1 79 1 E 

02 

-0  .2364755E 

02 

3 

0. 1238905E 

02 

0.1391628E 

02 

4 

0.4424074E 

01 

-0.7596874E 

00 

5 

-0.4208601E 

0 1 

—0. 28796  33  £ 

OO 

6 

-0. 2069723E 

00 

-0.  7998557E 

00 

7 

0 . 8464  32  7E 

00 

0.2104  794E 

01 

8 

0. 3291 799E 

00 

0.1254  503  E 

01 

9 

-0.557885AE 

00 

— 0.1  67  8884  E 

01 

10 

-0.2837964E 

00 

0.1566365E 

00 

I 1 

0.448E817E 

00 

0.  82591  ^2E 

00 

12 

-0 .6342903E- 

02 

-0. 1 633873E 

00 

13 

-0.  2921  963E 

00 

— 0.6565279E— 

01 

1 4 

0. 5386317E- 

01 

0.  1 49  9 777E 

00 

15 

—0.331 11 24  E— 

01 

0 .2785084E  — 01 

16 

-0.  1305802E 

00 

-0. 101 2938E 

00 

17 

0 . 547973 1 E— 

0 1 

0.54  1 2326E- 

01 

18 

0. 1251435E 

00 

0.191 8038  E 

00 

19 

-0. 208631 2E 

00 

-0.  1557431E 

00 

20 

0.1 179140E 

00 

0.1029936E 

00 

21 

-0.  1 1451  89E 

00 

-0.9353477E- 

01 

22 

0 .2414525E- 

03 

0.3330456E- 

Cl 

Note:  All  quantities  are  in  x,y  units. 
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0.3830773E-04 

65 

-0 . 1 16927 1 E—  0 2 

-0.  1 068  0 9 6E  — C3 

66 

0. 4905240  E—0  4 
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-0. 1 20  669EE  — 05 

1 05 

-0.83458948-06 

—0.5846961 E— 06 

106 

O. 31159648- 05 

0.  322  024  6 E — 05 

107 

-0. 2 1 67556  E—  0 5 

— 0.2048  7S6E— CS 

108 

-0.  1715321  E—  0 6 

-0.1  886759  E — 06 

109 

0. 9170A21E- 06 

0.  1 148447E-05 

110 

0. 1282769E-06 

0 .13475218-06 

111 

-0. 2130699E- 05 

— 0.  2 05 723 7E  — 05 

112 

0 • 13661 56E-05 

0.  1406525E- CS 

113 

0.  11289598-05 

0.12340  05E— 05 

114 

-0. 1028800E- 05 

- o.  i oi  aesiE-  05 
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I 


n 

B 

n 

C 

n 

I 15 

-0.131 344  85-05 

-0.1 246  324E-05 

116 

0.  1 5569635—  0 5 

0.  1 56R967E-05 

1 17 

-0 .B43634QE-08 

0.671 68S7E- 07 

1 IB 

-0.  752393  8C— 06 

-0.7546471 0-06 

1 19 

-0.  1622793E- 06 

-0.1  1 13204E-C6 

120 

0.1139069 F—  0 5 

0.1 15C875E-0S 

121 

-0.  6527602E-06 

-0.  61  68775E-06 

122 

-0 .7496735E-06 

-0.73927C5E- 06 

123 

0. 7258S87E- 06 

0.7580750E-06 

124 

0 • 522  1 764ET-  C6 

0.  5233  002E-  06 

125 

-0.961289SE-06 

-0  .9  29  4 666E-06 

126 

0. 200901  OF- 06 

0.2025572E  — 06 

127 

0 . 3 0992 5 8 E—  06 

0. 331664  IE- C6 

128 

-0.  106  762  5F—  06 

-0.1 01 7B22E-06 

129 

-0. 280262  OE— 06 

-0.264S621E-06 

1 30 

0. 2269661 E-06 

0.230  7490E-06 

131 

0. 9356683E- 07 

0. 1 0B5546E— 06 

1 32 

-0 . 1395750E-06 

-0.1 3H94 1 6E- C6 

133 

-0.  1 529258E—  06 

-0.1  394 1 74 E-06 

134 

0. 2 155629E- 06 

0.21681 C3E—  C6 

135 

-0. 34  0884  3 E— 07 

-0.2498253E-07 

136 

-0.  24 8794 4E-  0 7 

-0.227781 dE-07 

137 

— 0.7241 448  E— 07 

— 0.648258RE— 07 

1 38 

0.  6702049E-07 

0.6657  988  E — 07 

139 

0.  13881 2 IE—  0 7 

C.200e024E-07 

140 

—0 .41 86880  £—0  7 

-0.4101 790 E— 07 

141 

-0.  3901  64  1 E—  07 

-0.331 1 435E-07 

142 

0. 1 001300E-06 

0. 1 005498E- C6 

143 

-0. 4257138E-07 

—0 . 3833  354 E—  07 

144 

— 0. 4e05924E-07 

-0.471  6878E-07 

145 

0. 7 135057 E- 08 

0.1 083863E-07 

1 46 

0.  744  76  08 E—  07 

0.74937£3E— 07 

147 

— 0 • 53 74 5 7 BE—  07 

— 0.5108661E— 07 

148 

-0.2760153E— 07 

-0.2660B81E-07 

149 

0. 3 739 26 OE—  0 7 

0.  3956553E-07 

150 

0 . 1614653E-07 

0.1 6 2798  OE  — 07 

151 

-0.  S44902BC-07 

— 0.52  04  399  E— 07 

152 

0. 270 104BE- 07 

0. 2 7223  7 IE— 07 

153 

0.1 567296  F—0  7 

0.171500  7E  — 07 

154 

-0.  121 574  4E—  07 

-0.  1 1 68664E-C7 

155 

— 0 • 21 49 563 E—  07 

-0.201 B866E-07 

156 

0.  198  84  8 7E—  07 

' 0.2019971 E— 07 

157 

0. 3930207E- 08 

0.4938983E— 08 

158 

-0, 772594  4 E— 0 8 

-0.739A34eE-0fl 

159 

-0. 921 604 2 E—  08 

— 0.8195819E— 08 

160 

0.1 577 1 68E—  0 7 

0. 1577933E-07 
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APPFNDIX  C 


CONFORMAL  MAPPING  COEFFICIENTS 
FOR  THE 

EAST  COAST  REGION  FROM  PAMLICO  SOUND  TO  PENOBSCOT  BAY 


6 

= 0.2185243F 

02 

B 

o 

= C.  1 045535E 

03 

X 

= 374.0 

n 

B 

C 

n 

n 

1 

0. 84531 77E  01 

—0.7  3907 1 0 E 

02 

2 

-0.2573375E  02 

0.11 4 6656E 

02 

3 

0 • 2456  3 77  F 01 

—0.131 6729E 

02 

4 

— 0. 4 02  8893E  01 

0.80B9012E 

01 

5 

-0.1294892E  01 

— 0.580  3 839E 

01 

6 

0. 751 7704E  00 

0.2271 722  E 

01 

7 

-0.6833593E  0 C 

— 0. 606  82  85E 

00 

B 

0 . 87992  76  E 00 

0 .884  9853E 

00 

9 

-0.56B9548E  00 

-0.1988623E 

00 

10 

0.4  70985 2E  00 

0.31 53988E 

00 

11 

-0. 35051 02E  00 

—0 .140  5933  E 

00 

12 

-0.5908696F— 01 

-0.161  1297E 

00 

13 

0.1 1 06679E  00 

0.2047055E 

CO 

14 

-0. 3838463E— 0 1 

-0.7570612E- 

01 

15 

-0. 204S463E- 0 1 

0.3443124E- 

02 

16 

-0 .2993307  E— 0 1 

— 0.4593338E— 

01 

17 

0. 55191 73E- 01 

0.9201282E- 

01 

18 

-0 .3588944 E- 01 

-0.688^453E- 

01 

19 

0. 252  8971 E—  0 1 

0.5296R47E- 

01 

20 

-0. 2555390E- 0 1 

— 0.4  39624  6E  — 

01 

21 

0.21 58327F-0  1 

0.320  1 135E- 

01 

22 

-0. 2e8  36 1 6E—  0 1 

— 0.  31  50356E- 

01 

Note: 

All  quantities  are 

in  x,y  units. 
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V 


n 

B 

C 

n 

n 

23 

0. 1704723E-0 1 

0.19591  18E-01 

24 

-0.  821 84  82 E—  02 

-0.7724181 E— 02 

25 

0. 5179290E- OP 

0.30321  19E-02 

26 

-0  « 30  1 57  1 5 E-0  2 

0.1962418E— 03 

27 

0. 7938925E-02 

C.  4886530E-02 

28 

-0 .7356405E-02 

— 0.4150771E— 02 

29 

0.  212 03 03 E—  02 

0.1 651 054  E — 03 

30 

-0. 2383103E- 02 

-0.11  4 5 61 6E—  02 

31 

0.  120  64  59  E— 0 2 

0.4586522E-03 

32 

-0. 4140623F- 03 

0.  1 6660  64  E — 03 

33 

0 • 374  9528E— 0 3 

0.131 8797E- C 3 

34 

0. 4151922E-03 

0.5995615E— 03 

35 

-0. 6242518E- 02 

-0.6  71  6528E-03 

36 

0*4518004  E— 0 3 

0.4600859E- 03 

37 

-0.  1 1 79293E—  03 

— 0.8762796E— 04 

38 

—0* 1 105091 E—  0 2 

-0. 1 055562E- 03 

39 

—0*261 8267  E—0  3 

— 0.2268526E— 03 

40 

0. 33287 08E- 02 

0.  2 95 065  8E  — 03 

41 

-0 • 2192205  E—  0 3 

— 0.16891 96E—  03 

42 

0.  1223145E-  03 

0. 837  0787  E — 04 

43 

-0.820 1750  E— 04 

— 0.3827  059E— C4 

44 

0*  1276652E-0  3 

0.963981 3E-04 

45 

-0. 74  0 80  3 7E—  0 4 

— C. 4502  C9  0E—  04 

46 

0*4022771 E—  0 4 

0.207871  IE-04 

47 

-0. 65351 SOE—  04 

-0.51 0533OF-04 

48 

0- 6436008E-O5 

0.  1 1 4 9e4  0E— OS 

49 

—0 • 7533055E— 05 

-0.420  3695E-05 

50 

0. 2087624F- 04 

0.  3067982E-04 

51 

-0*  8448817E-05 

— 0.8  3950  30E— 05 

52 

-0.  2737144E-05 

— 0.2279463  E— 05 

53 

0. 109549ee-04 

0. 96294  6 3E  — 05 

54 

-0*6501 159E-05 

—0 .4  20 483 1 E — 05 

55 

-0.  5406307E-05 

-C.  72  332'*? OE  — C5 

56 

-0 .2243007E-05 

— 0.3883919E— 06 

57 

0.  423  8569E—  05 

0. 2996 593 E— 05 

58 

-0*  1688697E-  05 

- 0. 521 733 IE—  06 

S9 

0 • 35222  37  E—  05 

0.2865336E-05 

60 

-0*  2636244E-  05 

-0.2091  638E-05 

61 

0 • 146  808  7 C—  0 5 

0.  127 161 8E—  05 

62 

-0.  166440  IE-05 

— 0.1S47839E— 05 

63 

0.  1622032E-  05 

0.  1669C1  CE-OS 

64 

-0. 1230003E-05 

— 0.1246244E— 05 

65 

— 0. 2369671 E—  06 

-0.1 536769E-06 

66 

0*e677358E-07 

0.468 1 078E- 07 

67 

0. 7202404E— 06 

0.7946224E— 06 

68 

-0. 67964B3E-  06 

- 0.  7270P4PE-  C6 
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n 

B 

C 

n 

n 

69 

0 .43622596-06 

0 .49674636-06 

70 

-0.  1 076482E-06 

-0.1 4 08421 E — 06 

71 

0 .24590 1 9E- 0 7 

0.684  3044E- C7 

72 

-0. 83540276-07 

-0.10906096-06 

73 

-0. 6969583F- 07 

-0.40825536-07 

74 

0.87153536-09 

-0.1 3609706-07 

75 

-0.  1 0971  126-07 

0.32000266-08 

76 

0. 79048286- 0 7 

0. 75806216-07 

77 

0.26043766-0  7 

0.30559796-07 

78 

-0. 481 77236- 07 

— C.47AC5756  — 07 

79 

-0. 1 02.3679  6-0  7 

-0.95493856- 08 

80 

0.  3521  1 896-07 

0.37724276-07 

81 

-0.52718496- 07 

-0.54265596- 07 

82 

0.221  1 1 78F-07 

0.2568e4  7E—  07 

83 

-0. 76247206-08 

-0.95432466-08 

84 

0.  12521 96E-07 

0.15460926- 07 

85 

-0. 49960756-08 

-0.67233636-06 

86 

0 . 6 0990 3 6E- 0 8 

0.81 1 ie58E-08 

87 

-0.61197236-08 

-0  .72554026-08 

88 

0. 6175922E- 09 

0.  1 84 4 9 COE  — 08 

89 

-0. 3125592E-08 

-0.36191696-08 

90 

0.  560471  86-08 

0.6055163  E— 08 

91 

-0.61587006- 08 

-0.61235516-08 

92 

0. 1393683E-08 

0*13971 1 0 E—  08 

93 

0. 214  3876E— 06 

0.  23  8 84  7 6E  — 08 

94 

-0  .33712856-09 

-0.48583876- C9 

95 

-0.  20138346-09 

-0.35220486-10 

96 

0. 12193086- 05 

-0.5434  7606-10 

97 

-0.44785356-09 

-0  .2261 494E-09 

98 

-0.  26221  876-09 

— 0.36B781  5E  — 09 

99 

0 . 1 34  9 96  1 6-  09 

0. 27473226- C9 

1 00 

-0. 80127526-1 C 

-0.1 470834  E— 09 

101 

0.  14345386- 05 

0.21  1 S756E-09 

1 02 

-0.79446316-10 

-0.97774196-10 

103 

0.46998016-09 

0.49273986-09 

104 

-0 • 39365626-09 

-0.384  1 8656-09 

1 05 

-0. 68826366-1  0 

-0.71 951 72E-10 

106 

0. 29809146-  1 C 

0.47375236-1 0 

107 

0.68260866-1  0 

0 .58992236- 10 

1 08 

-0.  11512736-09 

— 0.9696301E  — 10 

109 

0.10110296-09 

0.90925716- 10 

1 10 

-0. 33274746-1  0 

-0.1 965071 E— 1 0 

1 1 1 

0. 24313456- 1 C 

0.  1 71  l 625E-1  0 

112 

-0.89472846-1 1 

-0.85250 16E- 12 

1 13 

-0.  5 0617616-12 

-O. 38813526-1 1 

l 14 

-0. 18969806- 1 0 

-0.14909556-10  ' 

I 


n 

B 

C 

n 

n 

1 15 

-0. 709 15  7 1 E— 1 1 

— 0 ,7886334 E — 1 1 

1 16 

0.  1821453E- 1 C 

0. 1 94  2 668E— 1 0 

1 1 7 

- 0 . 398  74  85  E— 1 1 

-0.31  1 5765E-  1 1 

118 

-0.  926062  3E-1  1 

-0.95321 16E-1 1 

1 19 

0 .7707409C-  1 1 

0.8761 047E-  1 1 

120 

0. 94  8 728  4 E— 1 1 

0.8745027E— 1 1 

121 

-0.  128127  2E— 1 C 

-0.1  1 531  68E-  1 0 

1 22 

0 • 54 17571 E— 1 l 

0.470  35 1 9E— 1 1 

123 

-0. 4854271 E-l 1 

— 0.387  0245E  — 1 1 

124 

0. 1 272995E- 1 1 

0.8C34556E- 12 

125 

— 0.8162463E— 1 2 

— 0.2376758E— 12 

1 26 

0.  2 4 38  756E—  1 1 

C.2222751E  — 11 

127 

-0. 1 364  84  7 E— 1 1 

— 0.107  8558E— 1 1 

1 28 

0.  714  6 806  E—  1 2 

0.6671 4 I5E-I2 

129 

-0.431 8938E- 1 2 

— 0.3278873E— 12 

1 30 

0. 1 00094  6 E— 1 1 

0. 1055302E-1 1 

1 31 

- 0. 20S2629E-  1 1 

-0.2062S67E-1  1 

1 32 

0.6150933  E— 1 2 

0.7 127636E- 12 

133 

0.  474  0078E-1 2 

0.4288399E— 12 

1 34 

-0. 393S754E- 1 2 

- 0. 3039050E- 12 

135 

0. 173488SE-1 2 

0.1 29 4 62 5 E— 1 2 

136 

0.  12266  1 3E- 1 2 

0.  1 895801  E — 1 2 

1 37 

-0 • 1 6249  64  E—  1 2 

-0. 1929000E- 12 

130 

0. 237601 8E-1 3 

0.6448195  E— 1 3 

139 

0.  1584382E-  1 2 

0.348401  OE  — 1 4 

140 

— 0 . 75754  29  E— 1 3 

— 0.5707148E— 13 

141 

-0. 1601686E-1 3 

— 0.1  69  8 1 61  E — 1 3 

142 

-0 .2932933E- 1 3 

-0. 2245293E- 13 

143 

0. 1-4569  76E-1 2 

0.1498009  E— 1 2 

144 

-0.  1 147427E-  1 2 

-0.  1 142878E-12 

145 

0. 29504  59  E— 1 3 

0,3  58  2 1 6 IE— 13 

146 

0. 2668761 E— 1 3 

0.2435500E-1 3 

147 

-0.2074897E- 1 3 

-0. 1477572E- 13 

148 

-0. 2806369E-1 3 

— 0.3055096E  — 13 

1 49 

0. 2675362E- 1 3 

0.  3 1 0 89 8 BE—  1 3 

150 

-0.247861 3E- 13 

— 0.2643  752E— 1 3 

151 

0.  1 65  0332E— 1 3 

0. 1933177E-1 3 

152 

-0. 122640BE-  1 4 

-0.21 20201E- 14 

1 53 

0. 7071394E— 1 4 

0.9649775E-14 

154 

-0. 9688173E- 1 4 

-C.99M  3S5E-14 

155 

-0. 349 79 9 4 E— 1 5 

0.37511 07E—  15 

156 

0. 411 24  3 0E—  1 4 

0.42S8614E-14 

157 

— 0 * 3 37625  OE—  1 4 

— 0.314346  9E  — 14 

158 

-0. 3620927E- 1 4 

-0. 33069 2 3E— 14 

159 

0. 5159257E- 1 4 

0.5141 R15E-14 

160 

-0. 6439003E- 1 5 

-0.31 03656E- 15 
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n 

B 

C 

n 

n 

161 

-0. 1550099E-1 4 

-0.1  644  377F.-14 

162 

0.  1889396E- 1 4 

0.21  76  993E-  14 

163 

-0.1462261 E-l 4 

—0 .1562678E— 14 

1 64 

0. 955 101 5F— 1 6 

0.  30851  69E-1  5 

165 

-0 . 186  361  BE-  1 5 

— 0.25635  39E-  15 

166 

0. 4924783E-1 5 

0 .52  32  3 75  E- 1 5 

167 

-C. 69751 1 7F- 1 5 

- 0.721 3222E- 15 

168 

0.20731 16E-15 

0.270  3232E  — 15 

169 

0. 2038455F— 1 5 

0.213991 7E  — 1 5 

170 

0 • 267  70  3 1 E- 15 

0.28  74  0 0 4E—  15 

171 

-0. 6432034E-1 5 

-0.61 62772E-15 

172 

0. 4125446E- 15 

0.4  083  4 02E  — 1 5 

1 73 

-0. 1303287F-15 

— 0.9925271E— 16 

174 

-0.  1 122934E-1 5 

— 0.1216983E— 15 

175 

0 . 4706550E- 16 

0,703437  6F  — 16 

1 76 

0. 55621 12E- 16 

0.4932 992E-16 

177 

-0. 728751 7F- 1 6 

— 0.5844  Cl  5E- 16 

178 

0,6455897  E— 1 6 

0.6  172  39  IE-  16 

179 

-0.  1810812  E—  l 7 

0.71 76773E-17 

180 

-0 , 251 9457E-  16 

-0. 2643164E- 16 

1 81 

-0. 258335  3F— l 6 

— 0.2060503E— 16 

182 

0. 9397427E- 1 7 

0. 9356655E- 1 7 

183 

0 • 284 99 6 5 E—  1 6 

0.31 53878E- 16 

1 e4 

-0. 435331 7E-  1 6 

-0.4277073E  — 16 

185 

0 .2342945F-  1 6 

0.2463311E— 16 

1 86 

0.6163340E-1 7 

0.774 3982E-1 7 

187 

-0.  1 038789E—  16 

— 0.  10571  59E  — 1 6 

188 

-0. 254  5895  E— 1 7 

— 0.5732719E— 18 

189 

0. e5822!6E-l 7 

0.8081 919E-17 

190 

-0. 1292874E- 16 

-0.11 26496E- 16 
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APPENDIX  D 


CONFORMAL  MAPPING  COEFFICIENTS 
FOR  THE 

GULF  COAST  REGION  FROM  LAGUNA  MADRE  TO  MARSH  ISLAND 


3 = 0*211 095  3E  02 

B = C.  9507664E  02 
O 

A = 360.0 


n 

B 

C 

n 

n 

1 

-0.4349387E  02 

-0.661 7682E  02 

2 

— 0*  .3 1 66200E  02 

0. 194  2 72  6E  02 

3 

0. 1242251 E 02 

— 0.1  17  382 1 E 01 

4 

-0.S063650E  01 

-0.1 447072E-01 

5 

0.1671185E  01 

—0.1921 633E  01 

6 

-0. 3597688E  0 1 

—0.673  3062  E 00 

7 

0.  1A19286E  0 1 

0.1366718E  Cl 

8 

— 0 • 1 179494  E 01 

-0.3091359E  00 

9 

0.4U9087E  00 

0.1  335779  E 00 

10 

-0.3775935E  OC 

0.4  832  629E—  02 

1 1 

0 *797 1597  E— 0 1 

0 .5321220 E—  0 1 

12 

-0.  1090030E  00 

0.1605034E  00 

13 

-0  * 1290084E-0 1 

-0.15598ee£  00 

14 

-0.  S 963 6 97 E—  0 1 

0,1 403255  E 00 

15 

-0.  1379531E-01 

- 0.63873  05E—  01 

16 

-0. 3508482E-0  1 

0.9618145E— 01 

17 

0. 66  04  981 E— 02 

-0.4127430E-01 

18 

-0 .14301  17E—0  1 

0.21 3 3995E—  0 1 

19 

0. 9 39  2757 E—  0 2 

— 0.2296Q28E— 01 

20 

-0.  1 132046E-  Cl 

0.2392460E— 01 

21 

-0. 1077687E-02 

-0 .2269067E— 02 

22 

-0.  1 27493 1 E—  0 1 

0.632  0797E-02 

Note:  All  quantities  are  in  x,y  units. 
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V 


n 

B 

C 

n 

n 

23 

0. 4 066765  E—  0 2 

0.141  1 9 75  E — 02 

24 

-0.  16674Q7E-  0 2 

0.4639465E- 03 

25 

0 • 399559  7 E— 0 3 

-0 .26  3 6 588E - 03 

26 

-0.  188621  9E—  0 2 

0.5682681  E — 03 

27 

0.9077224E— 03 

0.76  0 3890E— 03 

28 

-0. 1559198E-02 

-0.4942T92E-03 

29 

0.9190331E- 03 

0. 334  89C8E- 03 

30 

-0. 1 1 4 382  5 E— 0 2 

-0 .1035753E-03 

31 

0.  3 79523  0E-  0 3 

0.1  78  1 5 756.'  — 03 

32 

-0.67! 2531E-03 

0.948782  3E  — 04 

33 

0. 4371286E— 03 

—0.5  75  6 764  E — 05 

34 

-0. 583528eE- 03 

0. 1470665E- 03 

35 

0. 1 835 1 59  E—  0 3 

-0 .321 2969E- 04 

36 

-0.  2612983E-03 

0.304 1 208E-03 

37 

-0. 1043564E-04 

— 0.24  4 49 99E— 03 

38 

0.  102228  7 E-04 

0.232  69Q4  E — 03 

39 

-0. E843900E- 05 

-0.  1 970086E-C3 

40 

-0.  127881 3E-04 

0.1 54 1 862E-03 

41 

-0.  261  3426E-  04 

— 0.9325938E  — 04 

42 

-O. 14402958-04 

0. 103344 HE  — C3 

43 

—0.50798238— 05 

-0.5277937E-04 

44 

-0.  238344  7f  - f 

0.32466C9E-C4 

45 

0.26851098- 

-0.2672389E-05 

46 

-0.  431  1 55  1 r 

-0.2036  345E-05 

47 

0 . 26351 9 

0.  1 1 4P255C- 04 

48 

-0. 2369922? 

0.2308075E-05 

49 

0.  1 06  04  75E “ 04 

-0.281  666  OE—  05 

50 

-0.  U90996E-04 

0 .5396716E— 05 

51 

0. 79691 1 9E- 05 

— 0. 7 34  2 735 E — 06 

52 

-0.9435032F-05 

0. 305  690 IE— 05 

53 

0.  359 214  4 E—  0 5 

—0. 2 80  2860  E — 05 

54 

-0. 3163786E- 05 

0.  5202  766E-05 

55 

0 . 80697  3 3E— 0 6 

— 0.3  204  935E  — 05 

56 

-0. 31561 97E— 05 

0.2976863E-05 

57 

0 • 20821 34E—  05 

— 0.310024  5E  — 0 6 

58 

—0.27802 1 5E— 0 5 

0 .7  37  535 3E— 06 

59 

0. 1 130891E- 05 

-0.  71  6 86 70E  — C6 

60 

-0. 6122396E- CC 

0. 1 480043E- 05 

61 

-0. 18451 71 E—  06 

-0.1 119  7 1 OE  — 05 

62 

-0. 336542  OE—  06 

0.1101 696E-05 

63 

0.8938491 E—  0 7 

-0.454  240C.;-  06 

64 

-0. 301 2036E-06 

0.52354  4 7E—  06 

65 

-0. 3673406E- 08 

-0.  3582571  T-C6 

66 

-0 . 14  1 4085E-06 

0.42440588-06 

67 

0.  73883  65 E“ 0 7 

-0.1 289455E-06 

68 

-0. 3226247C- 06 

0.5 886  0 758 - 07 
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n 


n 

B 

C 

n 

n 

69 

0 . 27314  74  E—  0 6 

O. 100  2326E-06 

70 

-0.  2802S70E-  06 

— 0.3  232754E-07 

7 1 

0. 1234  1 09  E- 0 6 

0.1 08928CE- 07 

72 

-0. 949871 2F-07 

0.82R 1 825  E - 0 7 

73 

0.  1 907564E- 0 7 

— 0 • 6 229  69  ?E  — 0 7 

7 A 

-0. 3855900E-O  7 

0. 7623788E— 07 

75 

0. 1227572E- 07 

- 0.  4 72534  5E-07 

76 

— 0 • 2280 785  E—  0 7 

0 .595223 1 E— 07 

77 

-0.  4 7725  84  E — 0 0 

-0.391 7647C-07 

70 

-0. 1 30  59  3 IE-  0 7 

0.404  1 384E- 0 7 

79 

0 . 1 34  6654  E—  0 7 

-0. 1 214022E-07 

no 

-0. 21 78457E- 07 

-0.4  01  6400E-1  0 

8 1 

0. 28290  2 8 E- 0 7 

0.1231 37CE- 07 

82 

-0. 2833683E— 07 

— 0.554  2 695  E — 08 

83 

0.  13901 58E-  07 

0.4292090E- 08 

84 

-0. 1 041  305  E- 0 7 

0 .36289 3 4E  — 08 

85 

0. 2942560E- 0 e 

— 0.39004H2E-OB 

86 

-0  *2494349  E—  0 8 

0.7090  147E- 08 

07 

-0*  214  0656F-  08 

— 0.602374 OE— 08 

88 

0 . 1 858878E- 0 e 

0.  8 J9949CE- C8 

89 

— 0.3792454E— 08 

— 0.69248 18E— 08 

90 

0.  1 902798E- 08 

0.609681  9E-08 

9 1 

-0 .9630992C-09 

-0. 2973249E- C8 

92 

-0.  1 54  5032E—  0 8 

0.1 536916E-08 

93 

0,  17949ePE-  06 

0.3  73  0956E— C9 

94 

— 0 • 24545 63  E—  0 8 

— 0.555041 fcC—  09 

95 

0. 1989445E- 08 

0.9392  092E— 09 

96 

-0.201 8450E- 08 

— 0.632729  6E  — C9 

97 

0.  1354043E-08 

0 ,7495602 E— 09 

98 

-0.  1 189770E-  C 8 

— 0.  261  0C4  6E-09 

99 

0 . 560  70  7 7 E—  0 9 

0.922R99er- 10 

too 

-0.  2 8 726  8 OE—  09 

0.3  04  5 f\06F  —09 

101 

-0 • 989936  9E—  1 0 

— 0.3  76  64  85E— 09 

102 

0. 9587729E— 10 

0.5183951 E—  09 

103 

-0.  19201 87E- 09 

-0.3851  055E-C9 

104 

0 * 87 36652 E— 1 0 

0.3454421F—  09 

105 

-0.  652 1 29  9E- 1 0 

-0.21 78382E-09 

106 

-0.368493RE-  1 1 

0. 1826045E-C9 

1 07 

-0.31853 1 7E-  1 1 

— 0.8295703E— 10 

108 

-0. 58541 69E-  1 C 

0.  61 P8589F-1  0 

1 09 

0.5440021 E— 1 0 

-0 .706 1 680E- 1 1 

1 1 0 

-0. 7843175E- 1 0 

0.  3364A00E-1  1 

1 1 1 

0 .54853H9E- 1 0 

0.1 760203E-  10 

1 1 2 

-0. 5931984  E- 1 0 

-0.11  76  3 00  E— 1 1 

113 

0. 3508252E- 1 C 

0. 7755597C-  1 1 

1 14 

-0.3359514  E- 1 0 

0 .34384  77E— 1 1 
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- 


n 

B 

C 

n 

n 

1 15 

0.  1969359E-  1 C 

0.  369072  OF  — 1 2 

1 16 

-0 . 2062679F- 1 0 

0.5851 155F- 1 1 

l 1 7 

0.  1 077043E-1  0 

-0.161  1 202  E- 1 1 

1 10 

-0 . 1C97489E-  1 C 

0.677  09S4E-  1 1 

1 1 9 

0.3749065E-1 1 

— 0.5035324E— 1 1 

1 20 

-0.  3668546F- 1 1 

0.81 81 1 60E-1 1 

121 

-0 .4585235E- 1 2 

— 0.594  6987E-  1 1 

122 

-0. 7479945E-1 2 

0.74  30294  E-l 1 

1 23 

-0.  1297924E-  1 1 

-0.51 51255E- 1 l 

12a 

0. 93 3 1983 E- 1 3 

0 .535 1 474F- 1 1 

1 25 

-0.  851  2 74  5E  — 1 2 

-0. 34  6 569  6F  — 1 1 

1 26 

-0.462651  IE- 1 2 

0.3223064E- 1 1 

127 

0.  27 4 29  1 5F  — 1 2 

—0.1 365864 E-l 1 

1 28 

-0.  1541422E-  1 1 

0.8  7 8 2235E  — 1 2 

129 

0. 152029 1 E-l  1 

0.26  39  725E—  12 

130 

-0.  1 946227F-1  1 

— 0.28692P6E— 12 

1 31 

0. 1 42  878  OF-  1 1 

0.6296024E- 12 

1 32 

-0.  14  2 359 1 E— 1 1 

-0.1  924  273  E- 1 2 

133 

0. e095190E-  1 2 

0.24  1 622  IE- 12 

1 3a 

-0.6745650  E— 1 2 

0.1  149748E-  12 

135 

0. 2695209E-1 2 

— 0.1547703E— 12 

1 36 

-0 . 1687758E-  1 2 

0. 4 097  250E  — 12 

137 

-0. 992644  6 E— 1 3 

-0.3738514E-12 

1 30 

0. 690776  IE-  1 3 

C.4700165E— 12 

139 

-0. 1 29  4 38  3 E— 1 2 

-0.34060P1E- J2 

iao 

0.  1702639E-1 3 

0.291 1 647E-12 

la  1 

0.  1 289878E-  14 

-0. 1 232917E- 12 

1 42 

— 0.101 97  80E— 1 2 

0.90  33253E- 1 3 

143 

0. 9003627E- 1 2 

-0.6354  04  8E-14 

144 

-0.  1 2530t>9F—  1 2 

— 0,1 08  2 376E- 14 

1 45 

0. 1061812E-12 

0.4432039E— 13 

146 

-0.  11324 lOE-  1 2 

- 0.  1 36241  0E  - 1 3 

147 

0 . 70T  7363  E— 1 3 

0.2450861E— 13 

148 

-0. 6203187E- 1 3 

0.  3025605E-14 

149 

0 .3443273E- 1 2 

-0.1 103780E- 13 

150 

-0. 1978280C— 1 3 

0.37600P7E— 1 3 
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APPENDIX  E 


r 


•>< 

CONFORMAL  MAPPING  COEFFICIENTS 
FOR  THE 

GULF  COAST  REGION  FROM  MATAGORDA  BAY  TO  TIMBALIER  BAY 


B = 0.3266563E  02 

B = C.  64  0254  IE  02 
O 

X = 273.5 


n 

B 

c 

n 

n 

1 

-0.  1345455E  02 

— 0.1 26884  7E  02 

2 

— 0 . 9 308576E  01 

— 0.1253234E  02 

3 

0. 3063064E-0 1 

0.3363035E  01 

4 

0.  160121  3E  00 

0.1  09  1 622  E 01 

5 

— 0 • 58393  07E  OC 

0.3S96B8BC— 01 

6 

—0.1 9737 95 E 00 

— 0.2109445E  00 

7 

-0. 046  282  8E— 0 1 

0. 7170513E-02 

8 

0.546  74  1 IE- 0 2 

0.1 255086E  00 

9 

-0. 5171223E-01 

— 0.7066019E  — 02 

10 

-0.  1 0755O2F- 0 1 

-0.  110  0479E—  0 1 

1 1 

-0.980164 1 E-02 

0 .24728 C-2E— 01 

12 

-0. 6387778E- 02 

0.  ^ 62  8 1 63E  — 02 

13 

—0 .7022325 E—  0 2 

— 0 .2627  833E— 02 

14 

0.  14392 77E-02 

0.4072074  E-02 

15 

-0.  191 2029E-02 

0.15003  C7E—  02 

16 

-0. 1 5240  55  C—  0 2 

-0 .1732375E-03 

17 

0.  1 1 8022 1 E—  0 3 

0.  1 04  6393 E — 02 

18 

-0 . 1692594E-0  3 

0.3270761E— 03 

1 9 

-0. 24 28 7 62 E— 03 

0.1903720E— 03 

20 

-0. 3506050E- 03 

-0.1241541E-CA 

21 

0.10948256-0  3 

0.252445SE- 03 

22 

-0.  5571678E-04 

0.21 29799E-04 

Note : 

All  quantities  are 

in  x,y  units. 
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$ 

, 


n 

B 

C 

n 

n 

23 

-0  • 350  82 8b  E-  0 4 

0.404  04P8E- 04 

24 

-0.  3 78 1 4 96 E—  04 

0.1  14  4 81 6E-04 

25 

-0. 1 100774E- 04 

0. 7667904E— 05 

26 

0.114749SE-05 

0. 1734318E-04 

27 

-0.  1094296E- 04 

C.  1 076357E-05 

28 

-0.4207996  E- 0 5 

0.320  5 74  0E— 05 

29 

0.  51651  1 3E-06 

0.4804680E-05 

30 

-0. 4495374F- 05 

-0.  151  0 63  IE—  05 

31 

0. 1638793E-06 

0.20B0999E— 05 

32 

-0. 31 291 10E- 06 

0.  1 03  6 82  AE  - 05 

33 

-0 .639839 1 E-06 

0.  1503427C- 06 

34 

-0.  449821  5E-06 

0.241 6694E-07 

3S 

-0.  1 756886E- 0 7 

0.342421 8E—  06 

36 

-0.1 769853  E—  0 6 

0 .5766333E  — 07 

37 

0.  740571  0E-08 

0.  147  5 1 92E  — 06 

38 

-0 .756 1 528E—  0 7 

0.  16  7 3 0 2 7E— 07 

39 

-0. 774451 8E—  0 7 

-0.1 1 73998E— 07 

40 

0 • 698830 1 E- 06 

0.4872761E-07 

41 

-0.4358078 E— 0 8 

0.237O178E-07 

42 

-0.  1 5 03  878  E—  0 7 

0.263351 5E-08 

43 

-0. 1138452E-07 

0.4058880E- C9 

44 

-0. 2927226E— 08 

0.501 1 94 1 E— 08 

45 

-0.  1 103974E- Oe 

0.4151  303E-08 

46 

0 • 288  2630  E—0  9 

0.3492370E-08 

47 

-0.  292684  lE-Oe 

— 0.6412142E— 09 

48 

-0. 1887702E- 08 

-0.3581 659E- C9 

49 

0. 97331 3 9 E—0  9 

0. 1969457E-C8 

50 

-0. 3007952E— 09 

0. 34  54  06  8E  — 09 

51 

—0 • 79992 6 6 E— 09 

-0.361 3920E-09 

52 

0.  1161606E-09 

0.4086191E— 09 

53 

— 0 • 7 69 020 BE— 10 

0.  1 1 36903E-09 

54 

-0. 1 189104E-09 

0.428046eE- 1 1 

55 

-0. 2874996E- 1 0 

0.54361  14E-1  0 

56 

-0. 2802484E- 1 1 

0.547251 4E— 10 

57 

-0.  362  8853E— 1 C 

0.1S53833E— 12 

58 

-0.  1222868E- 1 C 

0. 121 2419E- io 

59 

-0. 1351430E- 10 

0.3202698E-11 

60 

0. 3004B4SF- 1 1 

0.  1421  820E-1  0 

61 

-0 .8372742E-  1 2 

0.631  296  IE  * 1 1 

62 

-0.  1 001  31  OE—  1 0 

— 0.5179364E— 11 

63 

-0.  154B850E-  1 1 

0.  1694  557E-  1 1 

64 

0 • 27621 7 4 E— 1 1 

0.4953277E-11 

65 

-0. 150501 1 E-l 1 

-0.6499526E-13 

66 

— 0 • 1 67704  Of  — 1 1 

— 0.7  356  5 6 HE— 12 

67 

-0. 515621 OE— 1 3 

0.6046279^—12 

60 

0. 67012B6E-  1 3 

0.4934  566E- 12 
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n 

B 

C 

n 

n 

69 

-0. 1041 3 35 E-l 2 

0 .1796432E-  12 

70 

-0.1 56Q783E- 1 2 

0.331 0604  E — 1 3 

71 

—0.1 83  75  09  E— 1 2 

— 0.50  1 4 89 8E— 13 

72 

0. 1567779E-1 3 

0.1 026  365  E — 1 2 

73 

0. 231SA09E-  1 2 

0.80538R7E— 13 

74 

— 0 . 5380629  E— 1 3 

-0.151 3999E-  13 

75 

-0. 203B72SE-1 3 

0.5883  84  7 E — 1 4 

76 

-0.320884  IE-  14 

0. 1 435 1 58E— 13 

77 

-0. 3887471 E-l 4 

0.7881 1 77E-14 

70 

-0. 371 1868E- 1 4 

0.401 5387E-14 

79 

—0 .271 00 63E— 1 4 

0.2644683E- 14 

00 

-0.  334 6254 E—  1 4 

0.2852854E— 15 

81 

-0 .2722021E- 15 

0.2089e6eE- 14 

82 

-0. 90542 19E— 1 5 

O.67S2058E-15 

83 

-0. 66864 36E- 1 £ 

0. 42  62  9 05E— 1 5 

84 

0 . 452  934  9 E— 1 6 

0.7726821E-  15 

85 

-0.  4 635597E—  1 5 

0.1653144E-16 

86 

— 0 • 4206397E— 1 5 

— 0.9398101E— 16 

87 

0. 1 195950E-1 5 

0.3406087E— IS 

88 

0. 2487004E- l 7 

0.  15297HE-15 

89 

-0. 1428360E-  15 

-0.431  0367E-  16 

90 

-0. 7437915E— 1 6 

— 0.7557572E  — 17 

91 

- 0. 9476944E- 1 7 

0.5468383E- 16 

92 

0 . 940  202 6E— 1 7 

0.40 18428E-16 

93 

-0. 1877212E- 16 

0.  1 51  06C2E-1  7 

94 

-0. 2430040E— 16 

— 0.1 058  724E— 16 

95 

0.  9067852E-1  8 

0.1 027978  E — 16 

96 

0.431517  0E- 1 7 

0.  1 06  0953E-  16 

97 

—0 .5278773  E— 1 7 

-0  .1 059649E—  17 

98 

-0. 3236399E- 1 7 

— 0. 3 837  736E  — 1 8 

99 

-0 . 1450895E— 1 C 

0.1789194E-  17 

100 

-0. 4036873E— 1 0 

0.B961303E— 18 

101 

-0. 4485G05E- 1 8 

0.4323934E— 18 

1 02 

-0.2379904  E— 1 8 

O .352  20  38E—  1 8 

103 

—0.370401 3E— 1 8 

O.311907OE  — 19 

104 

-0. 1238881E- 1 8 

0. 1467032E-  18 

1 05 

-0. 46 144 70 E-l 9 

0.1  364687E-18 

106 

-0. 5354953E- 19 

0. 6953575E— 1 9 

107 

—0,31 09555  E— 1 9 

0.5270040E- 19 

108 

-0. 5224460E— 1 9 

0.41  58  915E— 20 

109 

-0 .3378779E- 19 

0.444447  7E  — 2 0 

HO 

0.  1 095373E— 1 9 

0.3675055E— 19 

111 

-0. 64287 18F- 2 C 

0. 1 1 05266E— 1 9 

112 

-0. 151 7489E-19 

— 0.331 4 663E—  20 

113 

-0. 4 3 93 5 96 E—  2 0 

0.  35  85 689 E — 20 

114 

— 0.780607  3E- 2 1 

0.460 1 044E-20 
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n 

B 

C 

n 

n 

115 

-0.  1171447E-20 

0.2500088E-20 

1 16 

-0.7940360E-21 

0.  1 69324  9E  — 20 

1 17 

-0. 1 9 3352  0 E — 2 0 

-0 .261 0540F- 21 

118 

-0.  88  7 0977E-  2 1 

0.2520655E-21 

1 IP 

0 .2674543E- 2 1 

0. 1042381E-20 

1 20 

-0.  1754S23E-2  1 

0. 3496690E-P1 

121 

-0. 3950036E- 2 1 

— 0.  4 00  1 4 04 E-  22 

1 22 

-0. 166  2 322  E- 2 1 

0.7  3 9<}5  00E—  22 

123 

-0. S264034E— 22 

0.1 1 07630F-21 

1 2* 

-0.9163815E-23 

0. 1 029530E- 21 

125 

-0.2624163E-22 

0.4904399E— 22 

126 

-0.  6267675E-22 

-0.  1 1 73602E-22 

12? 

—0  *217  3399  E—  2 2 

0.1347  39 BE—  22 

128 

0. 3621 932E-23 

0.2753237E-22 

129 

-0. 741 7667E- 22 

0.8646540E-23 

1 30 

-0.8227958  F— 2 3 

0 .2775776E— 23 

1 31 

-0. 450  28  3 IE-  2 3 

0. 3098042E-23 

1 32 

-0 .4149869E-23 

0. 1 01 4040E- 23 

1 23 

0. 5361 73SF-2S 

0.3  56  7 867  E— 23 

134 

-0.2502105E-24 

0.21  4 5 265E—  23 

135 

-0. 18B18S3E-23 

-0 .2450 1 73E-24 

1 36 

-0.  95281  38E-24 

0.  1759  776E-24 

137 

0 .2794993E-2E 

0.790 1 839E- 24 

1 38 

-0.  126850  8E— 2 4 

0.3927325E— 24 

139 

-0. 2404606E- 24 

0.  1230278E-24 

140 

— 0 • 22574  53  E—  24 

0.30  2 552PE  — 25 

141 

-0. 9371 61 3E— 2 5 

0.7837279E-25 

142 

-0.61 1915  IE- 26 

0.11201 74E- 24 

1 43 

-0. 3691 736E-25 

0.4537704E— 25 

144 

-0. 4 0 7020  7E- 2 5 

0.7270641E-26 

145 

-0. 1 85  309  3 E- 2 5 

0.1 869  38 IE— 25 

146 

-0.  135821  IE-25 

0.126791  IE-25 

147 

-0.909264  5L- 26 

0.90772  16E-26 

148 

-0.  1953950E  — 26 

0.1 107780E-25 

149 

-0. 569P583E- 26 

0.31 05713E-26 

150 

—0 .7086726  E—  2 6 

-0.101 5449E-26 
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APPENDIX  F 


NUMERICAL  ANALOGS  OF  SURGE  EQUATIONS 

Admissible  finite  difference  analogs  of  equations  (40)  and  (41) 
which  are  nominally  centered  at  i , j , n (bearing  in  mind  the 
storage  of  variables  in  Figure  36),  are  as  follows: 

[Qs*(i> J >n+l)  - Qs*(i,j,n-lj/2At 

- f [QT* (i , j ,n+l ) + QT*(i,j ,n-l)]/2 

+ (gD/F(i, j)p(i))  [u(i  + l , j ,n)  - HB(i+l,j,n) 

- H(i-1 , j ,n)  ♦ HB(i-l,j,n)]/2AS* 

= xs*(i,j,n)  - (KqQC  i , j ,n-  1)/D2)QS#  (i , j ,n-*-l)  (F-l) 

and 

[QT*(i,j ,n*l)  - QT*(i, j,n-l)]/2At 

♦ f[Qs*(i,j,nfl)  + Qs*(i, j ,n-l)]/2 

♦ (gD/F (i , j)v(j)) [H(i, j+1 ,n)  - Hg(i , j+1 ,n) 

- H(i, j-1 ,n)  ♦ Hg(i. j-1 ,n)]/2AT* 

= tT*(i,j.n)  - (KQQ(i,  j .n-U/D2 )QT„ (i,  j ,n+l)  , (F-2) 

where  Q is  as  defined  by  equation  (54)  and  D as  defined  by 
equation  (61)  is  the  arithmetic  average  of  the  four  values  of  D 
about  the  point  i,j  at  which  the  flow  is  evaluated.  All  terms  are 
spatially  centered  at  i,j  and  all  but  the  bottom  friction  terms 
are  exactly  centered  at  time  level  n . The  latter  involve  flow  com- 
ponents at  the  new  level  n+1  and  at  the  old  level  n-1  within  the 
flow  magnitude  term  Q . This  form  is  known  to  lead  to  a stable 
algorithm  (Reid  and  Bodine,  1968).  The  Coriolis  terms  are  exactly 
centered  at  n as  well  as  i,j  . Equations  (F-l)  and  (F-2)  repre- 
sent two  equations  in  the  two  new  Q components  Qc*  (i , j ,n+l ) and 
Qj*(i,j,n+1)  in  terms  of  quantities  at  previous  times  and  are 
readily  solved  for  each  component  individually.  The  resulting 
explicit  relations  for  these  components  are  given  by  equation  (55) 
where  Gj  , G2  , G3  are  defined  in  equations  (58)  to  (60). 

The  centered  finite  difference  analog  of  the  continuity  equa- 
tion (42)  leads  directly  to  the  explicit  relation  for  H at  the  new 
time  level  given  by  equation  (57). 
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APPENDIX  G 


MODEL  VERIFICATION 

Prior  to  the  surge  computations  the  conservative  properties  of 
the  numerical  model  were  examined.  The  model  was  applied  first  to 
seiches  in  an  enclosed  square  basin  12  kilometers  on  a side  and  5 
meters  deep.  The  number  of  grid  increments  along  each  side  was  20. 
To  invoke  the  most  severe  seiching  mode,  the  water  surface  was 
initially  assigned  a linear  slope  downward  from  one  end  of  the  basin 
(j=.JM)  to  the  opposite  end  (j=l).  The  initial  setup  for  the  nonro- 
tating basin  computations  was  50  centimeters,  but  this  was  reduced 
to  10  centimeters  when  rotation  was  allowed  to  suppress  unnecessary 
fluctuations.  The  time  step  taken  for  all  test  runs  was  50  seconds. 

The  simulation  of  a uninodal  seiche  in  a nonrotating  basin  was 
considered  first.  With  the  specified  time  step,  the  period  of  the 
severest  mode  is  theoretically  68.57  At. 

Figure  G-l  shows  that  the  computed  fundamental  period  is  within 
one-time  increment  of  the  theoretical  period. 
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Figure  G-l.  Computed  water  level  at  grid  point  (21,20) 
for  uninodal  seiche  test  case. 


The  time  variation  of  the  volume,  potential,  kinetic,  and  total 
energies  are  given  in  Figure  G-2.  Each  of  these  quantities  are  area- 
weighted  sums.  This  evaluation  scheme  should  exhibit  small  perturba- 
tions more  readily  than  unweighted  sums  of  the  variables  over  the  grid 
The  variation  of  the  volume,  potential,  kinetic,  and  total  energies  is 
small,  with  the  potential  energy  showing  the  largest  deviation. 
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Rotational  effects  on  the  seiching  mode  were  checked  by  specify- 
ing the  Coriolis  parameter  as  0.75  o0,  1.00  oQ,  and  1.25  a0  where  aQ 
is  the  angular-  frequency  of  the  fundamental  mode.  The  computed  fre- 
quencies, o , were  compared  with  those  obtained  by  Rao  (1966)  and 
Platzman  (1972). 

Table  G-l  gives  the  comparison  with  Rao's  and  Platzman's  results. 


Table  G-l.  Comparison  of  computed  values  of  o/oQ  with  those  obtained 
by  Rao  (1966)  and  Platzman  (1972). 


f/°o 

O 

D 

D 

Average 

computed 

o/°o 

% error 

From  Rao 
(1966) 

From  Platzman 
(1972) 

0.75 

0.769 

0.76 

1.1 

1.00 

0.723 

0.721 

0.75 

4.0 

1.25 

0.686 

0.683 

0.71 

3.9 

Although  the  errors  occurring  with  f/aQ  > 1.00  are  large,  these 
values  of  f/oQ  represent  stern  tests  of  the  algorithm.  It  is 
definitely  encouraging  that  no  evidence  of  instability  is  indicated 
in  the  simulated  hydrographs. 
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APPENDIX  li 


WIND  DEFORMATION  PROCEDURE 


Assuming  that  the  wind-stress  components  (T$*  , ty*  ) 

are  known  in  the  stretched  shelf  coordinate  system  from  the  applica- 
tion of  a symmetric  wind  field  model,  the  problem  in  part  is  to 
identify  a region  for  a given  hurricane  and  coastline  where  it  is 
appropriate  to  alter  these  stresses  to  reflect  the  influence  of  land. 
The  proposed  deformation  equations  are  applied  only  at  points  which 
are  located  within  this  region.  Existing  charts  of  hurricane  winds 
and  the  investigation  by  Graham  and  Nunn  (1959)  provide  the  basis  for 
the  empirical  deformation  formulas.  In  this  manner,  the  analytical 
representation  obviates  the  detailed  input  of  a massive  sequence  of 
digitized  wind  field  data  that  conforms  with  the  observations  near 
the  coast  as  well  as  offshore. 

T,  . * . „ (SYM)  , (SYM) 

The  wind-stress  components,  tg*  and  Xy*  , are 

altered  at  a point  for  which  Yj  > Yp  according  to  the  relationship: 


(SYM) 


(H-2) 


where  m is  a constant  chosen  to  be  2 and  1^  is  an  influence 
factor  given  by: 


if  - 1 - d£  r (Yj  - y/v2  * (h-3) 

where  R^  is  a constant  taken  as  the  average  radius  to  maximum  winds 
of  the  hurricane,  Df  is  a distortion  factor,  Yp  is  the  shortest 
distance  the  point  is  from  land,  and  Yj  is  the  distance  the 
influence  region  extends  from  the  coast  relative  to  the  point  in 
question  (Figure  H-l).  In  general,  Yj  should  depend  on  Xp  , the 
distance  the  point  in  question  is  from  a centrally  located  point  along 
the  coast  (xc,yc).  This  latter  point  is  generally  assigned  the 
coordinates  of  intersection  of  the  hurricane  track  with  the  coastline. 
From  Figure  H-l,  the  distances  Xp  and  Yp  are  given  by: 

XD  = ±(  (xs  - xc)2  + (ys  - yc)2]!*  , (H-4) 

and 


* 


(H-5) 


vD  -l(Xp  - *s>!  • <yp  - y,)2l'‘  , 

where  (xs,yg)  is  a point  along  the  coast  determined  by  the  inter- 
section of  a line  normal  to  the  coast  and  passing  through  the  point 
in  question  (*p>yp)-  The  sign  convention  is  that  Xp  is  positive 
if  the  point  in  question  is  to  the  right  of  the  hurricane  track  and 
negative  if  it  lies  on  the  left.  In  general,  the  influence  region 
is  defined  by: 

YI  = [ W + Sf  1 ' Rf  ’ (H"6) 

where  Fj  is  the  influence  function  of  Xp  , Sf  is  a shift  factor 
which  permits  the  influence  line  to  be  moved  normal  to  the  coast, 
and  Rf  is  a range  factor.  The  influence  function  attempts  to  re- 
late the  degree  of  distortion  in  the  wind  isovels  relative  to  the 
tract  of  the  hurricane  and  the  shape  of  the  coastline.  The  extent 
of  the  land  influence  is  given  by  Yj  , which  depends  on  Sf  and 
Rf  . These  coefficients  as  well  as  Df  are  considered  to  be 
related  to  the  distance  of  the  hurricane  from  the  point  where  it 
crosses  the  coast. 

Additional  investigation  is  required  to  establish  a generalized 
procedure  based  on  physical  principles  for  deforming  the  winds  in  a 
manner  which  faithfully  reproduces  the  observations.  Consequently, 
each  storm  is  parameterized  individually  to  accomplish  the  task.  The 
relationships  used  for  deforming  the  winds  of  Hurricane  Carla  are: 


FI  * 1 XD  - \ l/?  ’ 

(H-7) 

Df  - .1  f 1 ♦ •-•05<Dh/V2>2]  , 

(H-8) 

Sf  - M .5  - if.-'  (Dh/2RhM  ^ . 

(H-9) 

r 1.5  , if  D,  < -SR, 

n h 

• 

Rf 

. .5  . e-025'VVS)2- if  Dh  i$5h 

(H-10) 

where  R,  is  taken 
n 

as  approximately  37  kilometers  and 

Dh  ■ * 1 <*c  - V2  * O'c  - V21'1  • 

(H- 11) 

The  sign  convention 
water  and  negative 

is  that  is  positive  if  the  hurricane 

if  the  storm  is  over  land.  For  Hurricane 

is  over 
Camille, 

v 
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dl 

_ — : — 

f 

these  relations  are: 


(XD  ♦ 6RJ/I.S  , if  XD  < -31^ 

| XD  |/1.5  , if  -3^  < XD  < 0 ( H - 1 2 ) 

| XD  |/7  , if  XD  > 0 

Df  * 0.2  , CH- 13) 

Sf  = 0.5  , CH- 14) 

and 

Rf  = 1.0  . (H-15) 

The  Hurricane  Gracie  wind  fields  were  not  deformed  because  the 
winds  are  presented  in  digital  form.  Any  deformation  obtained  would 
be  dependent  upon  a subjective  analysis  which  requires  additional 
developmental  work.  Moreover,  the  selection  of  Gracie  was  made  on 
the  knowledge  that  the  results  would  possibly  be  compared  with  those 
obtained  by  other  investigators  using  symmetric  wind  fields. 


APPENDIX  I 


SYMBOLS  AND  DEFINITIONS 
Transformation  coefficients 

Nondimcnsional  wind-drag  coefficient 

Depth  of  water 

Distortion  factor  for  deforming  the  wind 

Distance  of  hurricane  from  the  coast 

Maximum  expected  depth  of  water  at  any  point 
within  computing  grid 

Depth  of  water  relative  to  mean  sea  level 

Error  function 
Coriolis  parameter 

Scale  factor  associated  with  the  orthogonal 
curvilinear  coordinate  system 

Influence  function  for  deforming  the  wind 

Acceleration  due  to  gravity 

Water  level  relative  to  mean  sea  level 

Hydrostatic  elevation  of  the  sea  surface  corres- 
ponding to  the  departure  of  the  atmospheric 
pressure  from  a constant 

Influence  factor  for  deforming  the  wind 

Number  of  computational  grid  points  along  the  S*  and 
T*  axes 

Wave  number  = ir/X 

Wave  number  of  free  wave  in  the  annulus 
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rh 
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S*  ,T* 
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V„ 


wc,ws 


»10 


Nondimensional  bottom-drag  coefficient 

Number  of  Fourier-type  transformation  coefficients 
in  the  series 

Atmospheric  pressure 

Far  field  pressure 

Central  pressure  of  the  hurricane 

Volume  transport  per  unit  width  in  the  S*  and 
T*  directions 

Polar  coordinates 

Distance  from  the  hurricane  center  to  any  point 

Range  factor  for  deforming  the  wind 

Distance  from  the  hurricane  center  to  the  region  of 
maximum  winds 

Shift  factor  for  deforming  the  wind 

Distance  normal  to  the  seaward  boundary  and  along 
the  coast 

Coordinates  of  the  stretched  shelf  coordinate 
system 

Long  wave  travel  time 
Maximum  winds  of  the  hurricane 

Weighting  factors  for  the  coast  and  seaward 
boundary  curves 

Windspeed  at  an  elevation  of  10  meters  above  the 
water  surface 

Rectilinear  coordinates 
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Vyc 


vyh 


xp’yp 


x ,y 
s s 


c c s s 
X ,Y  ,X  ,Y 


Z 

0 

A 

At 

£ 


c s 
Y » Y 

X 


♦ 

PA’PW 


°c*  >°n 


Coordinates  of  landfall  of  the  hurricane 

Position  of  hurricane  center 

A point  in  the  influence  region 

A point  along  the  coast  determined  by  the  inter- 
section of  a line  normal  to  the  coast  and  pass- 
ing through  (xp,yp) 

Coordinates  of  the  coast  and  seaward  boundary 
curves  to  be  mapped 

Extent  of  the  land  influence  for  deforming  the 
wind 

Plane  of  (x,y) 

Extent  of  n (±)  in  the  ?-plane 
Arc  length 

Surge  algorithm  time  increment 

The  alongshore  curvilinear  coordinate  of  the  shelf 
coordinate  system 

Prorating  factors  for  the  length  of  the  transform- 
generated curves 

Half  the  horizontal  extent  of  the  region  to  be  mapped 
in  the  Z or  ?-plane 

Functions  transforming  £ and  n to  the  stretched 
shelf  coordinate  system  (S*,T*) 

Ingress  angle 

Density  of  air  and  water 

Total  and  nth  contribution  to  the  "curvi linearity" 
variance 

Bottom  resistance  stress  divided  by  in  the  S* 

and  T*  directions 
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Wind  stress  divided  by  pw  in  the  S*  and  T* 
directions 

Angle  between  the  £ and  x axes 

The  offshore  curvilinear  coordinate  of  the  shelf 
coordinate  system  (£,n) 

Plane  of  U,n) 
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